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Abstract 

We investigate tlie use of noncommutative Grobner bases in solving partially prescribed 
matrix inverse completion problems. The types of problems considered here are similar to those 



in [BJLW95]. There the authors gave necessary and sufficient conditions for the solution of a 
two by two block matrix completion problem. Our approach is quite different from theirs and 
relies on symbolic computer algebra. 

Here we describe a general method by which all block matrix completion problems of this 
type may be analyzed if sufficient computational power is available. We also demonstrate our 
method with an analysis of all three by three block matrix inverse completion problems with 
eleven blocks known and seven unknown. We discover that the solutions to all such problems 
are of a relatively simple form. 

We then perform a more detailed analysis of a particular problem from the 31,824 three 
by three block matrix completion problems with eleven blocks known and seven unknown. A 



solution to this problem of the form derived in [BJLW95] is presented 



Not only do we give a proof of our detailed result, but we describe the strategy used in 
discovering our theorem and proof, since it is somewhat unusual for these types of problems. 



'University of California, San Diego, Math Dept. 9500 Oilman Dr., San Diego, CA 92093-0112. Partially supported 
by the AFOSR and the NSF. 
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1 The problem 



We consider block matrix completion problems similar to that in [BJLW95|. Here, we take two 
partially prescribed, square matrices, A and B, and describe conditions which make it possible to 
complete the matrices so that they are inverses of each other. That is, we wish the completed 
matrices to satisfy 

(1) AB = I and BA = I. 
1.1 A sample problem 

An example of such a problem is: given matrices fci, ^2, ^3, and ^4, and 

(2) A=[ "2 \ andS=(' 

is it possible to find matrices ui, U2, M3, and M4 such that equation (|^) is satisfied? The answer to 



this question, due to [BJLW9f;], is given in Section 2.3 



We now describe our problem in detail. 

1.2 The general block matrix inverse completion problem 

We begin by partitioning two matrices, A and B, whose entries are elements in an arbitrary field F, 
conformally for matrix multiplication into nhy n block matrices. Next, we choose / of these blocks 
to be known and 2n^ — / to be unknown. We give some conditions on the known matrices, which may 
be expressed algebraically, such as invertibility or self-adjointness. We will now define our problem. 

We ask if it is possible to fill in the 2ti? — I unknown blocks so that equation (0) is satisfied and 
seek to derive formulas for these matrices in terms of the prescribed blocks. To be more specific, 
we might even call this problem the purely algebraic partially prescribed matrix inverse completion 
problem. The solution to such a problem will be a set of matrix equations in the known and unknown 
submatrices. 

2 The solution 

In general, it is not known how to solve a system of matrix equations where several of the matrices 
are unknown. Unknown matrices can appear in matrix equations in several ways, of which some 
are more computationally acceptable than others. We will analyze these forms and classify certain 
solution sets. In this section, as well as throughout the paper, noncommutative variables ki will be 
considered initially known, and noncommutative variables Ui will be considered initially unknown. 
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We first recall some standard definitions. If F[A;i, ... , fc;, ui, . . . , U2n2-;] is a noncommutative, 
free algebra over the field F of characteristic 0, then a subset 2 will be called an ideal if /, g e X 
implies that If frf + Iggrg £ X for lf,rf,lg, and Vg arbitrary elements of F[fci, . . . ,ki,ui,... , U2n2-/]- 
Given a finite set of polynomials {pi,... ,Pm} in F[/si,... ,fc;,ui,... ,U2n2-;], we say that X is 
generated by {pi, . . . if X is the smallest ideal containing {pi, . . . ,Pm}- 

2.1 A good, triangular solution 

The following definition is useful in identifying a class of problems, particularly those consisting 
of matrix equations, which are usually computationally tractable. This will be demonstrated in 
Section 3.5.2| . It is, in general, impossible to verify condition (||) in Definition 2T below. For this 



reason, we give two versions of our definition: a weaker, computable, version and a stronger, in 
general incomputable, version. We will later introduce approximations to condition which can 
be verified. 

Definition 2.1 Let X be an ideal in ¥[ki, . . . ,ki,ui, . . . ,U2n2-/]- We say that X can be made 
weakly formally backsolvable if there exists a bijective map 

cr :{!,... ,2^2 ^ {!,... ,271^ - /} 

and a finite set G of polynomials which generates X such that 

G = GoUGiU...UG2„2_, 

where Go — G C\ F[fci, ... , /c/] and 

Gi C F[A;i, . . . , ki,Ucr(i),- ■ ■ ,u„{i)\ \F[fci, . . . , A:/, u<^(i) , . . . ,u<T(i-i)], G,; ^ 0, 

for 1 < i < 2n^ — L //, in addition, we have the condition, 

(3) no proper subset of G generates X, 

we say X can be made formally backsolvable. We say that the set G is formally backsolvable. 



Definition 2.1 , for many people, will be more intuitive in an expanded notation. What follows is an 



intuitive and expanded notation for weakly formally backsolvable. Indeed, the set G of polynomials 
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J Q2n^~i.m„ 2_,} t^^*^ form 

(4) qo,i{ki,... ,ki) = 

(5) go,2(fci, • ■ ■ , fci) = 

(6) go,mo(fci, • ■ ■ = 

(7) ,fc;,u^(i)) = 

(8) ,fc;,u<^(i)) = 

(9) g2,i(fci,--- ,fci,u^(i),u<^(2)) = 

(10) q2,m2{kl, ■ ■ ■ ■.kl,Ua{l),U„{2)) = 

g3,l(fci, . . . ,fci,U„(l),U^(2),U^(3)) = 

g2n2_i_i,m^^2_^_^(fci, . . . ,A:;,Uo.(i),Uo.(2), ■ ■ ■ ,Ucr(2n2-l-l)) = 
92n2-i,l(fcl: ■ • ■ , fc;, Ua-(i), Uo-(2), Uo-(3), W(t(4), ■ • ■ 7 Uo.(2n2-l)) = 

(11) 92n2-i ) = 

where the ki are known, rrii > 0, the are unknown, and cr is a permutation map on integers 1 
to 2n^ — I. The first subscript of the polynomials gr,™ indicates the number of unknowns allowed in 
the equation, as well as implying the existence of the bolded unknown Uo.(r)- 

We refer to equations (|[^) which contain only knowns as compatibility conditions on the 
knowns. These equations, in only the known variables, must hold if a completion is possible. 
Equations ([7[]ll|) containing unknowns we call equations triangular in the unknowns due to 
the triangular structure exhibited by the unknown variables. 



Usefulness for computation is discussed in Section 3.5.1 
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2.2 Decoupled solutions 

The following definitions are useful in identifying two classes of problems, particularly those consist- 
ing of matrix equations, which are usually even more computationally tractable than the formally 
backsolvable form, introduced in the previous section. We call the two classes essentially decoupled 
and formally decoupled. It is not always possible to verify condition (jlj) in Definition 2^ below. 



For this reason, we give two versions of each class definition: a weaker, computable, version and 
a stronger, in general incomputable, version. Later, we will introduce approximations to condition 
(|l|), which can be verified. 

Definition 2.2 Let X he an ideal in ¥[ki, ... , fc;, wi, . . . , U2n2-i] • Let there exist j , I < j < Ir? — I; 
an injective map a : {1, . . . , j} {1, . . . , 2n^ — I}; a bijective map r : {j + 1, . . . , 2n^ ^ ^ 
{1, . . . , 2n'^ ~l}\image{a) ; a set G* of polynomials; and a finite set G of polynomials which generates 
I such that 

G Go UGi U . . . UGj U {ur[i) - gr{ki, . . . ,kuu„i^i), . . . ,m<,(j)) : j + 1 < i < 2n'^ - 1} 

UG* 

where Gq = G H F[fci, ... , fc;] and 

G, C F[fci, . . . , fc;, \ F[fci, . . . , fc;], G, ^ 0, 

for 1 < i < J. 

We now make the following definitions concerning the ideal introduced above. 

1. We say that the polynomial ideal I can be weakly essentially decoupled and that the set 
of polynomials G is weakly essentially decoupled. 

2. If G* = 0, then we say that the polynomial ideal T can be weakly formally decoupled and 
that the set of polynomials G is weakly formally decoupled. 

An important nondegeneracy condition is 

(12) no proper subset of G generates X. 



3. If condition ( 12 ) holds, then we say that the polynomial ideal X can be essentially decoupled 



and that the set of polynomials G is essentially decoupled. 



If G* = and condition (12) holds, then we say that the polynomial ideal X can be formally 



decoupled and that the set of polynomials G is formally decoupled. 
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Definition |2.2| , for many people, will be more intuitive in an expanded notation. What follows is 
an intuitive and expanded notation for weakly decoupled. Indeed, the set G of polynomials, 

G — {90, 1-90,2, • • ■ , 9j,mj , + l - Ur(j + 1), ■ ■ • 1 <ls j +2n'^ -l-j " ^ir(2n2_;)} 

in the formally decoupled case or 



G — {go.l, 90,2, ■ ■ • , 9j,mj , qsj + l - Mr(j + 1), • ■ ■ , 1sj+27i^-l-j - ""r (2^2 -/) , 9sj +2^^ + 1 , • ■ ■ , ^sq } 

in the essentially decoupled case, has the form 

(13) qo.iih,... ,ki) = 

(14) qo,2{ki,... ,ki) = 

(15) <?0,mo(^l; ■■■ ,ki) = 

(16) ,fc;,U^(i)) = 

(17) ,fc;,UCT(i)) = 

(18) g2,i(fci,.-. ,fci,u^(2)) = 

(19) ,A:z,u^(j_i)) = 

(20) . . . ,fc,,u^(j)) = 

(21) . . . ,fc,,u^(j)) = 

(22) "r(j+l) = 9sj+l(fcl,--- 

(23) U^(2n2-l) = 9sj+2n2-;-j(fcl, ■ ■ • , fc;, . . . ,Mcr(j)) 



T/ie following equations, which we call compatibility conditions on the unknowns, 

will not exist in the formally decoupled or weakly formally decoupled cases, but might 
occur in the essentially decoupled case. 
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(24) 



9s3+2n2-/-j + l(fci, . . . , A;;, Ml, . . . , U2n2_i) = 



(25) qs„{ki,... ,ki,ui,... ,U2n^-i) = 

In the above system of equations sj — X]i=o with rrii > and sq —\G\. The first non-zero 
subscript of the doubly subscripted polynomials g^,™ indicates the occurrence of one and only one 
unknown Uo.(r)- 

Equations of the form (|2^p3|) will be referred to as singletons. A singleton equation is character- 
ized by the fact that there is a single instance of an unknown variable which does not occur in equa- 
tions (p^pT|) . This unknown variable appears in the singleton equation as a monomial consisting of 
only itself. The singleton variable is the left hand side of equations ( p2[]23| ), UT{j+i), ■ ■ ■ ,itr(2n2-i)- 
Singleton equations in Definition 2.2 are Mr(i) ~ 9i{ki, • ■ ■ , ^z, Mcr(i), ■ • • , Wo-q)) = 0. 



The singleton equation has a very attractive form for a human who wishes to find polynomials 
in few unknowns. Given an equation in knowns and unknowns £, it allows one to eliminate the 
unknown singleton variable, for example from £ by replacing instances of the unknown 

indeterminate with its equivalent polynomial representation, in the example case Qsj+i- After this 
substitution has been performed, the equation £ will not contain the unknown singleton variable. 

As in the formally backsolvable case, we have compatibility conditions on the knovi^ns 
(P^15). These equations, in only the known variables, must hold if a completion is possible. All 



unknown variables ui, . . . ,U2„2_; in equations ( [13[ - |25| ), which are not singleton unknowns, appear 
in equations ( p6[ - pl| ) without any other unknown variables. Therefore, we think of this system of 
equations as decoupled. We call equations (p^pT|) equations in one unknown. 

In the formally decoupled case, the coupling compatibility conditions on the unknowns (p^^5|) 
are absent. This is obviously a better form of solution than the essentially decoupled form, since 
any solutions for Mo-(i), • ■ • , Mo-(j) will do. In the essentially decoupled case, one must verify potential 
solutions for ■ • • T^cr{j) with equations ( p4[]25| ). 

Notice that these decoupled solution forms, essentially decoupled and formally decoupled, satisfy 
the formally backsolvable criteria. We have the following set inclusion relationship: 
Formally decoupled C Essentially decoupled C Formally backsolvable 
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2.3 A sample answer 



Here, we give the solution to the problem presented in Section |l . l| . In |BJLW95 |, it was shown 



that, for invertible fc^, the matrices A and B defined in satisfy (|l|) if and only if the unknown 
submatrix U4 satisfies the following relation 

(26) Uik2Ui = U4 + kskik^. 

The other unknown submatrices are then given in terms of M4: 

(27) Ml = k^'^ - k2Uik^^ 

(28) M2 = k^"^ — k'^^u^k^ 

(29) M3 = k/Ji-^u^k^^k^^ 

This answer contains no compatibility conditions on the knowns. Equation ( p6| ) is an equation in one 
unknown M4. The remaining equations (|2^p9|) are singletons. Therefore, the equations associated 
with this matrix completion problem can be formally decoupled. In such circumstances, we would 
say that the problem can be formally decoupled. In the language of Definition Go = 0, G* = 0, 
J = 1, a(l) = 1, Gi = {(pel)}, t(2) = 2, r(3) = 3, and r(4) = 4. 



This main theorem of [BJLW95 is simpler, from a computational perspective, than the results 



we are presenting here, since (0) contains fewer equations and fewer variables. In addition to being 



proven via traditional methods, the main theorem in |BJLW95 was also proven using noncommu- 



tative Grobner methods in |HS9!;]. 



3 Main results on 3x3 matrix inverse completion problems 

We have performed extensive analysis of the 3x3 block matrix inverse completion problem. In 



particular, we have concentrated on the problem described in Section 1.2 where n is three and I is 
eleven. We have assumed in our detailed analysis that all eleven known blocks are invertible. 

We begin by noticing that if one matrix problem is a permutation of another, then a solution 
to one transforms to a solution to the other. We then define a property that characterizes certain 
matrix completion problems which we call strongly undetermined. We will present a classification 
result which characterizes the solutions to our seven unknown, eleven known matrix completion 
problems. In this result, strongly undetermined problems are the worst behaved class. 

This section also includes a detailed result on a particular matrix completion problem, which is 



in the same spirit as the | BJLW95 result described in Section 2^. We will conclude this section 



by outlining a method by which one may find a numerical solution to a matrix completion problem 
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using our symbolic solution. Much of the work in this paper appears in the Ph. D. thesis of Dell 
Kronewitter, [KroOO . 



3.1 Configurations and permutations 

In our investigations of 3x3 block matrix completion problems, we will refer to a configuration as 
a classification of blocks into knowns and unknowns. We will specify a configuration with k's and 
It's. For example, 



is a configuration. 

A 3x3 permutation matrix is a 3x3 matrix consisting of three I's and six O's. No two I's may 
appear in the same row or the same column. There are, of course, six such matrices. 

For a given 3x3 configuration of knowns and unknowns, one may apply 3x3 (block) permutation 
matrices, 11 and to A and B to get n^^^^I' and '^^^BH and obtain at most 36 other equivalent 
configurations. That is. 



if and only if 

U-^A^ ^-^BU = / and ^-^^BH n^^A* = /. 

In describing solutions A and B to this problem, we will only give one member of a particular 
equivalence class {Il^^A^,'i'^^BIl}. 

3.2 Strongly undetermined 

Assume that the pair of block matrices, A and B, are partitioned into known and unknown blocks 
that are compatible for matrix multiplication. 

Definition 3.1 A and B are said to be strongly undetermined if there exists an entry of the 
block matrices AB or BA, which is a polynomial consisting entirely of unknown blocks. 

Notice that A and B being strongly undetermined is equivalent to the existence of both an entire 
row (column) of unknown blocks in A and an entire column (row) of unknown blocks in B. For 
example, the following configuration of known and unknown blocks is strongly undetermined. 



(30) 




AB = I and BA = I 



(31) 
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The product of these two matrices, AB, has the following form. 

uifcy + M2fc9 + uafcio uiks + U2U5 + uskii U1U4 + U2U6 + U3U7 
kikr + k2kg + kskiQ kiks + k2U5 + kskn kiU4 + k2Ue + ksuj 
kj^k-j + k^kg + fcefcio k4^kg, + k^u^ + k^kn k^u^ + k^UQ + fcguy 

Since the upper right entry (in boldface) is a polynomial made up entirely of unknown blocks, 
configuration (pll) is strongly undetermined. 



3.3 A class of 31,824 3x3 matrix inverse completion problems 

In our investigations, we have analyzed (via computer) a certain collection of 3x3 matrix completion 
problems. Two 3x3 block matrices have a total of 18 entries. We have analyzed those which have 
seven unknown and 1 1 known blocks and do not have the strongly undetermined property. We have 
chosen to put efforts into this ratio of known to unknown blocks because we believe Theorem |2| the 
initial subject of our research, to be surprising, and yet lack the computational resources to study all 



3x3 matrix completion problems, or even one 4x4 matrix completion problem. Section 4.4 describes 
how the motivated researcher with unlimited computational power can go about analyzing a block 
matrix problem of any size of the type addressed in this paper. 

The following theorem shows that all of our seven unknown, 11 known block matrix completion 
problems (which are not strongly undetermined) have particularly nice solutions. 

Theorem 1 Let A and B he three by three block matrices such that 11 of the 18 blocks are known 
and seven are unknown. Let the known blocks be invertible. The corresponding partially prescribed 
inverse completion problems may be classified as follows. 

1. If the configuration of unknown blocks is not strongly undetermined, and is not of the form 
given in (pi) or a permutation of such configuration, then the partially prescribed inverse 



completion problem can be weakly essentially decoupled, in the sense of Definition 2.L 



2. Problem (|5£) can be made weakly formally backsolvable in the sense of Definition 2.1. Thus, 
all but the strongly undetermined cases can be made weakly formally backsolvable. 

These answers, that is the resulting weakly formally backsolvable or weakly essentially decoupled 
systems of equations, satisfy a technical non-redundancy condition, compatibility 3-nondegeneracy, 



which will be defined in Section J^.c once we have built up our Grobner machinery. 



The exceptional case mentioned in Theorem |i| is 

ki k2 ks \ / kj kg U4 

(32) A = I A;4 k^ kg \ and -B = kg kio 

ui U2 Us I \ fell ue wj 
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The proof of this theorem, which requires noncommutative symbohc software, wiU be given in 



Section 4.f; . Answers to the individual problems, which consist of sets of polynomials similar to that 
found in equations (|2^p9[), can be found on the internet at 
http://arXiv.org/abs/math.LA/0101245. 

3.4 Detailed analysis of a particular 3x3 matrix inverse completion prob- 
lem 

We now give a closer analysis than that given in the last section of a particular matrix inverse 
completion problem, which satisfies the assumptions of Theorem |l|. We show how someone, interested 
in a particular matrix completion problem, might arrive at a finer analysis of the problem instead 
of the rather terse conclusion given in Theorem |l|. Our goal in this section is to present a short, 
computationally simple set of formulas which give the solution to a particular partially prescribed 
inverse matrix completion problem. Our conclusions will have the same flavor as those presented in 



BJLW95I 





We will analyze a particular problem from those addressed in Theorem |l|, the known/unknown 
configuration: 



(33) A = \ u c w and B = 



or the equivalent permuted form. 



A 



where a through k are known and invertible block matrices, and the underlined t through z are 
unknown block matrices. 





b 


t \ 
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(= 


e 


w 


and S = 


J 


z 






V 


c / 




u 


i 


y J 



Theorem 2 Given A and B as in (|ff,j[) with invertible knowns a,b,c,d,e,f ,g,h,i,j , and k, as well as 
the invertibility of the matrix made up of the outer known blocks of A in 

a b 

e 



> d 
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then AB = I and BA = I if and only if the knowns satisfy the following compatibility conditions: 

(35) p{da-^-eb-^) = {a-^b-d~^e)q 

(36) {da-^ - eb^^y^q ^ (da-^ - eb-^)-^qe{da-^ - eb^^y^q 

+ {da^^ — eb^^y^qdg + jb(da^^ — eb^^y^q — kci + jag 

(37) p{a-^b~d-^e)-^ = p{a-^b - d-^e)~^ep{a-^b - d'^e)-^ 

+p{a^^b — d^^ey^dg + jbp{a^^b — d~^e)~^ — kci + jag 

where 

(38) p = {-a-^h-^i + a-^bjh'^i~d-^ -~d-^ejh-^) 

(39) q = {-kf-^a-^ + kf-^gda-^ -b-^ - kf-^geb-^) 

The unknown matrices can then be given as: 

(40) z = {a-^b - d~y)-^{-a-^h-^i + a-^bjh-^i ~ d'^ - d-^ejh'^) = (a^^fe - d-^^y^p 
or equivalently 

(41) z = {-kf-^a-^ + kf-^gda-^ - b-^ - kf'^ geb-^){da-^ - eb-^)-^ = q{da-^ - eb-'^)-^ 
and then 

(42) t = -agi^'^ - bzi~'^ 

(43) u = -k^^ja - k^'^zd 

(44) V = k-^ ~ k-\jb - k-^e 



ezi ^ 



(45) w = i^^ - dgi~^ 

(46) X = + fk-^j - gda-^ + fk-^zda-^ 

(47) y = c-^k-^jaf + c'^^k~^jbk + c-'^k'^zdf + c-^k-^zek 



This answer consists of three compatibihty conditions on the knowns (|35|-|37D , an equation in one 
unknown ( pO[ ) or (pj), and singletons (|4^^7|), and is, therefore, formally decoupled. 



The proof of this theorem wiU be given in Section 4.6. We mention some key points about how 
the proofs in Theorem |l| and Theorem |^ compare in Section |3.6| . Solutions to all of the problems 
(configurations) addressed in Theorem |^ can be found via the internet at 

http://arXiv.org/abs/math.LA/0101245. These solutions consist of a formatted list of equations 
in both I^lJ^pC and Mathematica form. 
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3.5 Numerical solutions to matrix equations 

Now we want to see how symbolic solutions to matrix completion problems may be applied to 
numerical problems in order to find numerical matrix completions. Let us see what is involved 
in the numerical solution of a matrix completion problem with configuration ( p3[ ) which was solved 
symbolically in Theorem 2. Assume, for example, that the matrices A and B are 12x12 and therefore 
each block is a 4x4 matrix. That is, we are given 4x4 matrices a, 6, c, d, e, /, g, h, i, j, and k and are 
trying to find 4x4 matrices t,u,v,w,x,y, and z which will form the completed inverse. We now 
apply Theorem 2 and the formally decoupled set of equations ( ^5[ ] 

The first step is to determine whether the matrices a, 5, c, d, e, /, g, h, i, j, and k are compatible. 
That is, they must satisfy equations (|3^p7|). If these equations are satisfied, then the next step is 
to determine a value for our unknown z using either equation ( ^0| ) or (|4l]). Since we have assumed 
the invertibility of a^^b — rf^^e, the coefHcient of z, this step cannot fail. Finally, one can determine 
values for the 4x4 matrices t,u,v,w,x, and y from the singleton equations (|4^p7|). If all of these 
steps have occurred successfully, then we have formed our inverse matrix completion. This illustrates 
very general behavior which we now describe. 

3.5.1 Solving a decoupled system of matrix equations 

A general decoupled system of matrix equations is made up of compatibility conditions on the knowns, 
equations ( ITsfjr^ ); equations in one unknown, equations (16-2l|); singletons, equations ( p2[p3| ); and 



possibly compatibility conditions on the unknowns, equations (|2J-^5|). 

Given such a decoupled set of matrix equations, one can first verify that a completion may be 
possible by verifying equations (p^p^ containing only the given (known) matrices. Then one can 
use equations ( p6[]2l| ) 

(714 = 0, . . . , Qj^rnj = 

to simultaneously solve for the (possibly non-unique) matrices ■ ■ ■ , u^-q) or to determine that 

solutions do not exist. Notice that this may constitute a difficult numerical problem by itself, 
especially if the matrices under consideration are of large dimension. It is then a simple matter to 
find matrices u^-(^j_^_l^, . . . , u^(2n^-\) by evaluating polynomials 

The boldface u's in equations ( |l6[]2l| ) indicate the unknown matrix being solved for in each step. 

In the essentially decoupled case, one must check that the solutions ui, . . . , U2„2_/ which this pro- 
cedure derives are acceptable by validating compatibility equations (p^]25|). Of course, an advantage 
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of decoupled equations is that their sohition may easily be parallelized. 
3.5.2 Solving a formally backsolvable system of matrix equations 

A general formally backsolvable system of matrix equations is made up of compatibility conditions 
on the knowns, equations (H^), and equations triangular in the unknowns, equations ([^p]). One 
can first verify that a completion may be possible by verifying the compatibility conditions on the 
knowns (||). 

Next, one attempts to use the equations triangular in the unknowns to solve for the Uk matrices. 
One may solve for the (possibly non-unique) Uo.(i) using equations (^[^) or determine that a solution 
for does not exist. Notice that this may constitute a difficult numerical problem by itself, 

especially if the matrices under consideration are of large dimension. With the results obtained for 
Ug.(i), one may next use equations (p[jlO|) to solve for Uo-(2) or to determine that a solution does not 
exist. This process continues until we have solved for all unknowns, that is until we have formed 
an inverse completion, or have determined that a completion is not possible. The boldface u's in 
equations (|7|-|ll|) indicate the unknown matrix being solved for in each step, if a solution can be 
found for each boldface matrix. 

3.6 The strength and Hmitations of our method 

The difference in strength between Theorem ^ which was derived automatically by computer algebra, 

and Theorem which concentrated on one case and employed some human intervention, illustrates 

the limitations of our automatic methods. Theorem reduced the particular completion problem 

it addressed, configuration (p3[), to solving a set of compatibility conditions on the knowns and a 

set of singletons defining each of the unknowns in the problem. This is of course the most desirable 

form of solution. On the other hand, Theorem |^ reduced all but 36 of the 31,824 matrix completion 

problems it addressed to essentially decoupled equations, a highly informative but less desirable 
answer. 

The way the stronger answer in Theorem ^ was derived illustrates the role of human intervention. 
First, we apply Theorem]^ to problem ( [33| ) and obtain an essentially decoupled set of equations E. 
One equation has the form 

{a~^b — dr^e)z = q 

where g is a polynomial. We assumed that a^^b—d~^e is invertible, and implemented the assumption 
by adjoining the equations defining the inverse to E and rerunning the Grobner basis algorithm.^ 
Naturally, this solved for z, in other words, produced a singleton defining z, and thus derived those 

^ Also there was a bit of beautification of formulas which was not essential to the form of the result. 
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compatibility equations on knowns resulting from substituting for z in the equations in E where it 
appeared. 

The only human intervention behind Theorem ^ was the decision that a^^h — dT^e is invertible. 
However, a human would not know that it is critical before Theorem ^ was applied. Without making 
such an invertibility assumption, the results in Theorem |l] are as far as one can go. 

In fact, the invertibility assumptions and subsequent computer manipulation, which transform 
equations (|6^|6^ into a singleton ( ^0| ) or (|4^) in z, are typical of the sort of human intervention 
which is required in many problems. 

4 Solving the purely algebraic inverse matrix completion prob- 
lem 

In this section, we will describe a method for solving general matrix completion problems of the type 
described above. The main tool we will use for our solution of the problem is the production of a 
noncommutative Grobner basis. We will review Grobner basis definitions and results, and present a 
pure algebra interpretation of our matrix completion problem. This section also contains the formal 
proofs of the results presented in Section |[ 

4.1 Background on Grobner bases 

Grobner bases are a useful tool in the manipulation and analysis of polynomial ideals. We will 
review how the Grobner basis may be used to 

I Discover whether a polynomial p is a member of a polynomial ideal /. 

II Show two polynomial ideals / and J are the same, given generating polynomial sets gi and gj. 

Ill Transform a set of equations into an equivalent set with a "triangular" form, described in 
Section ^. 

4.1.1 Monomial orders 

Essential to the polynomial machinery we use is the existence of a total order on the monic monomials 
in the polynomial algebra under consideration. 

Recall the definitions of lexicographic and graded (length) lexicographic monomial orders on 
commutative monomials, as discussed in 1CLS92f . The noncommutative versions are essentially 



similar, but to ensure a well defined total order a monomial mayQbe parsed from left to right in the 
tie breaking length lexicographic order criteria. 



^In fact, this is the scheme used in the NCGB computations. 
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The NCGB software uses a combination of these two types of orders, which we will find useful. 
It lets one define sequential subsets of indeterminates such that each subset is ordered with graded 
lexicographic ordering within the subset, but in which indeterminates of a higher set are lexico- 
graphically higher than indeterminates of a lower set. That is, a monomial consisting of one element 
in a higher set will sit higher in the monomial ordering than a monomial consisting of the product 
of any number of elements in lower sets. The NCGB notation uses the ^ symbol to discern the 
subset breakpoints discussed above. For example, when we write xi < X2 < ^ x^ we get that 
X1X2 < X2X1, X3X2X1 < X4, and X3 < x\X2- We call such an ordering multigraded lexicographic. 

4.1.2 Lead terms and Grobner rules 

Given a polynomial p, there exists a unique term of p whose monomial is highest in such an order. 
Denote this LeadTerm(p). For example, if we have X3 < X2 < xi then 

LeadTerm(xi — X2X3 + x^) — x^. 



For technical reasons, | CLS92 |, our orders must satisfy the following relation for any polynomials 
p and q, 

LeadTerm(p) LeadTerm((7) = LeadTerm(p • q). 

With this definition of leading term, we can introduce replacement rules. 

Every polynomial p corresponds to a replacement rule T(p), where the left side of the rule 
(LHS^RHS) is the leading term of the polynomial, and the right side is the negative sum of the 
remaining terms in the polynomial. If we have xi > X2 > x^, then the polynomial xi — X2X3 + x\ 
corresponds to the rule x\ X2X3 — xi. We may write 

r(xi - X2X3 + xl) = xl X2X3 - xi. 
The next example illustrates how we can apply a set of replacement rules to a polynomial. 

Example 4.1 To the polynomial 

a;ia:ia;2a;i + X2X3X1 + X2, 

we may apply 

{xl X2X3 - Xl, X2X3XI 4:Xl} 

to get 

{X2X3 - Xi)x2Xi + (4X1) -\- X2 — X2X3X2X1 - X1X2X1 + 4xi + X2- 
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A polynomial p may be reduced to a hopefully simpler form with a set of polynomials P by 
applying the replacement rules r(P). Reducing a polynomial, by a set of rules for commutative one 
variable polynomials, is similar to the classical Euclidean division algorithm one might use to find 

say ^ ■ scenario of this paper, it is easier to view reduction as a replacement scheme 

where rules are applied to a polynomial to change it to a "simpler" form. 

Example 4.2 The polynomial a::2Xia:2a;i + X2 may be reduced by the polynomial X1X2 + 3: 
The Euclidean division algorithm approach: 

(a;2a;ia;22;i + X2) - X2{xiX2 + 3)xi = -3x2Xi + X2 
The replacement rule approach: 

To X2Xia;22;i + X2 apply r(a;ia;2 + 3) = a;ia;2 —3 
to get — 3a;2a;i + X2- 

If repeated application of these types of rules to a polynomial transforms the equation to 0, 
then we have shown that the polynomial under consideration is an element of the (two-sided) ideal 
generated by the relations used to create the rules. 

4.1.3 Grobner bases 

A Grobner basis G for a polynomial ideal / enjoys the powerful property that a polynomial q is 
an element of the ideal if and only if repeated application of all rules T(g) arising from Q send the 
polynomial q to 0. See [CLS92 for the commutative version or | HS99| for the noncommutative 



generalization. One says that the Grobner basis solves the ideal membership problem, which is 
Problem I given in the beginning of this section. One may write 

Given generating sets 7W and TV, one can use this property of Grobner bases to show that the 
ideals (M) and (JV) are the same. Given Grobner bases Gm and GjV for M and A/", respectively, we 
will have 

{M) = (A/-) 



if 



m ^ for all polynomials rn E Ai and n ^ for all polynomials n G J\f. 
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Thus, Problem II, presented above can be solved, if finite Grobner bases can be found for both of 
the ideals. It may be solved, without computing a Grobner basis, by using the reduction properties 
of a partial Grobner basis, as will be done in Section 4.7.3| . 



4.1.4 Generating a Grobner basis 

If the indeterminates commute, then a Grobner basis is always a finite set of polynomials, and there 
exists an algorithm called Buchberger's algorithm which finds this set, given a generating set of 
relations for the ideal. In the noncommuting case, a Grobner basis for an ideal may be infinite. 



Nevertheless, there exists a similar algorithm due to F. Mora [Mor86|, which recursively defines a 
Grobner basis and terminates if it happens to be finite. In practice, even this finite Grobner basis 
may be incomputable when computer resources are taken into account. One stops the algorithm 
after a specified number of iterations, thereby generating some finite approximation to a Grobner 
basis. 

This finite approximation, though not exhibiting the omniscient powers of a true Grobner basis, 
is often useful in reduction, as will be shown below in Sections 4.7.2 and 4.7.3| . These finite approx- 



imations are, of course, made up of elements of the ideal generated by the original relations. There 



are a few available software implementations of this noncommutative algorithm (|HS97|, NCGB and 
| KG98[ |, OPAL). For our computations, we used the package NCGB. 



The form of this (partial) Grobner basis is dependent on the order under which Mora's algorithm 
is executed. One places variables which one wishes eliminated high in the order. The order 

fcl < fc2 < ■ • • < fc; ^ Ul ^ U2 ^ • • ■ <C ti2„2_; 

will cause the output of the Grobner basis algorithm to have the form of the equations ffllll) as much 



as is possible, as will be discussed in Section 4.1.5. We may therefore be able to solve Problem III 



introduced on page p7|, if a good enough approximation to a Grobner basis is practically computable. 
4.1.5 Elimination theory 

As promised above, we now motivate the multigraded lexicographic ordering with a central concept 
in elimination theory. 

Definition 4.3 A monomial order on the monic monomials in ¥[xi, . . . is said to be of the 

j-th elimination order if a monomial containing one of {xj+i, . . . ,Xn} is ordered higher than any 
monomial made up of the indeterminates {xi, ... ,Xj}. 

A Grobner basis Q for an ideal X, created under a j-th elimination order, exhibits the following 
ideal relation (C/P|F[a;i, • ■ • ,Xj]) = If]¥[xi, . . . ,Xj]. In words, this says that all polynomials re- 



20 



suiting from our original generating relations which contain only the variables {.ti, ... ,Xj} can be 
found in the ideal generated by C/p|F[a;i, • ■ • ,Xj]. This property, well known in the commutative 



case, was extended recently to noncommutative algebras in [HS99|, Theorem 11.3. We will use this 



property to assist with the triangular goal described in Sections 2.1 and 2.2 



4.1.6 A small basis 

Often, what a human wishes to find is not a set of relations with the reduction properties described 



in Section 4.1.3, but a small set of relations which describes the solution set. We refer to such a small 
set as a small basis. Such a goal may be accomplished by iteratively including the first k elements 
of our basis _B in a set Bk, creating a partial Grobner basis Gbk from Bk, and trying to reduce the 
other polynomials B\Bk with this partial Grobner basis Qb^ ■ When Bk has the property that the 
excluded relations B\Bk are elements of the ideal generated by the included relations Bk, our goal 
has been achieved. All of our relations lie in the polynomial ideal generated by Bk- We call such an 
algorithm the small basis algorithm. 

The sequence in which relations are presented to the small basis algorithm is obviously important. 
The small basis algorithm acting on (^x^ , x^ , x, l) returns the unenlightening (x^, x^, x, l) , but when 
presented with {l,x,x'^,x^) the algorithm returns (1). (Computer time required discourages some 
idealized implementation, which would consider all permutations of our relations.) 

4.2 A pure algebra interpretation of the purely algebraic partially pre- 
scribed inverse matrix completion problem 

This section describes our matrix completion problem in the language of an algebraist. The reader 
may skip this section, if so desired, with no loss of continuity to the paper. 

Labeling the known blocks as ki, we may consider the free algebra F[fci, ... ,ki] over the field 
F under consideration, modulo some presupposed conditions (e.g. the invertibility of a known 
submatrix ki, which is expressed in ideal theoretic notation as {kik^^ — 1, k^^ki — 1)) 

F[fci,... ,ki] 



(48) S : 



conditions on the knowns ) 



Let the size of the square submatrices be m. Picking the known block matrices consists of defining 
a map 

Defining 

(49) T = S[U1, . . . ,U2n2_i], 
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completing the matrices A and B may be viewed as a map 

$ : T ^ Af„,(F) 

such that 

$|5= (j). 

Here, we are interested in a special completion $ so that our matrices satisfy (|l|). Let Flatten be the 
operation which takes a set of matrices to their constituent blocks. If J is then the ideal generated 
by relations (j^), 

(50) J = {Flatten{AB -I,BA~ /)), 

our goal is achieved when J lies in the kernel of <i>. 
Given G, a generating set for J, so that 

(G) = J, 

our compatibility conditions on the knowns have been defined to be G n S. These are equations 
in the backsolvable case or equations (13-1^) in the decoupled case. Notice that all of the relations 



making up G n ( J \ S*) will contain at least one Ui. 
4.3 Nondegenerate solutions 

On closer analysis of the weakly formally backsolvable form described in equations (^-pl|) , the exis- 
tence of qo,i implies the existence of qi^i, ■ ■ ■ , q2n'^-i,m^ 2 , ■ Simply multiply go,i by the appropriate 
Uj. A more interesting set of relations has the following property. 

Definition 4.4 A set of equations {pj = : 1 < < s} will be called nondegenerate if 
(51) ^ 



Condition (^) in Definition 2.1 or condition in Definition 2.2 is equivalent to the nondegenerate 



property. Due to the infiniteness of the noncommutative Grobner basis, condition ( pl| ) cannot in 
general be verified. A condition which can be verified computationally is the following. 

Definition 4.5 A set of equations {pj = : 1 < < s} will he called f-nondegenerate if 

(52) P^ ^ {{PJ}J^^)^ 

where is the i iteration partial Grobner basis created from {pj}j^i. 
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Consider the system of equations 

(53) fci + fcafci = 

(54) fciMi + k2kiui = 

(55) kik2 + kiU2k2U2 + kiU2 

(56) k2k2U2 + Uik2ki = 

where ki and ^2 are known and ui and U2 are unknown. 

These relations appear to be essentially decoupled, but are only weakly essentially decoupled, 
since equation ( |5^ ) is a member of the ideal generated by equation (|5^). By removing equation 
(Q) we are left with a formally backsolvable system of equations, since we have an equation in 
one unknown U2 and an equation in two unknowns, ui and U2. Nondegeneracy can be verified by 
attempting, and failing, to reduce to each polynomial in |(|5^), (|5^ , (^6|)| using a Grobner basis 
created from the other two polynomials. For this example, the three Grobner bases are finite, so 
this procedure is definitive. 

4.3.1 The special case of compatibility nondegenerate solutions 

An approximation to nondegeneracy, often used in this paper, is compatibility nondegeneracy. This 
form of nondegeneracy is the condition that the equations which contain unknowns cannot be reduced 
to by those equations which contain only knowns. 



Definition 4.6 We will call a weakly decoupled set of equations of the form compatibility 
nondegenerate if 

(57) QhiM ^ {{qoA, • ■ • , 90, mo}) for hi > and qh ^ ({<7o,i, • ■ • , 90,™,,}) for h> Sj + 2n^ - I - j. 



That is, the relations (|16|-P5D which define the for i — 1, . . . , j are not trivial, and are merely 
consequences of the compatibility conditions on the knowns, equations (13 1^). The singleton equa- 



tions ( p2[p3| ), those which define Wr(i) for i = j + 1, . . . , 2n^ — I, are obviously not trivial, since the 
term Ur(i) cannot be reduced by any Grobner rule containing only ki, . . . ,ki. 
We also have the computational analogue. 

Definition 4.7 We will call a weakly decoupled set of equations in the form of ( jZ^ - p^ compati- 
bility ^-nondegenerate if 

(58) qhiM ^ ({90,1, • • • ■qo,mo})e for hi > and qh ^ ({90,1, • ■ ■ ,9o,mo})« for h> Sj + 2n^ - I- j- 
where ({(70,1 j ■ ■ ■ , <lQ,mo})i is the I iteration partial Grobner basis created from {qo,i, ■ ■ ■ , 9o,mo}- 
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A weakly formally backsolvable system of equations can also be compatibility ^-nondegenerate. 

Definition 4.8 We will call a weakly formally backsolvable set of equations in the form of 
compatibility f-nondegenerate if 

(59) qhiM ^ iUos, ■ ■ ■ ,qQ,m„})e for hi > 

where {{qo,i, • ■ • , 9o,mo})^ is the £ iteration partial Grobner basis created from {qo^i, ■ . ■ , 9o,mo}- 

Beware that these definitions are algorithm dependent, since the Grobner basis algorithm allows 
for some variability in how it is implemented. Furthermore, research into different variants of the 
noncommutative Grobner basis algorithm has not thoroughly addressed the reduction properties of 
partial Grobner bases. If one were to run the Grobner basis algorithm for an infinite number of 
iterations (which might result in infinitely many polynomials), then one could verify condition (^7|). 

Our computational resources are, of course, finite and we can do no such thing. Still, our three 
iteration partial Grobner basis offers a computational approximation to the condition (|57|). This 



form of non-redundancy given in Definition 4.7 was used to verify compatibility 3-nondegeneracy 
for all the problems in Theorem |l|, which were essentially decoupled. All but 36 of them were of 
this form. While this is all that wc did automatically on all 31,824 cases, Theorem ^ serves to show 
what one can do by further applying our Grobner basis methods to a particular case. In that case, 
we gave a concise solution to the matrix completion problem without an "infinite computation" . 

4.4 A recipe for solving the general block matrix inverse completion prob- 
lem 

We are given matrices A and B partitioned conformally for matrix multiplication into blocks 
each and a configuration of I prescribed (known) and 2n^ — I unknown blocks. We may also be given 
conditions on these matrices which are expressed algebraically (e.g. invertibility, aa^^ — 1 = and 
a^^a — 1 = 0). We look to discover compatibility conditions on the known matrices and formulas 
for the unknown matrices to solve our problem, that is, to ensure (^ is satisfied. 

This paper shows that this goal may often be achieved by following the steps below. 

I Fill in the known blocks of A and B with symbolic, noncommuting indetcrminates, fci, . . . ,ki. 
II Fill in the unknown blocks of A and B with symbolic, noncommuting indetcrminates, 

Ui, . . . ,U2„2-i- 

III Create the noncommutative polynomials resulting from the operations AB — / and BA — I. 
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IV Create a (noncommutative, partial) Grobner basis for the polynomials derived in step III and 
any assumed algebraic conditions on the matrices under the order: 

fcl < /C2 < ■ • • < fc; < Ul -C U2 < • • ■ -C M2n2_/ 



V Check that the result has some attractive form, such as those described in Section (|1^^5|) 
or Section l^l, (Hl|). 



VI Verify that the relations defining unknown matrices are not merely consequences of the other 
relations by using the Small Basis Algorithm or some variant of it. 

The noncommutative algorithms we use are not yet well understood and, therefore, their effectiveness 
on a particular class of problems can only be determined by experimentation. 

4.5 Proof of seven unknown, 11 known theorem 

We created a Mathematica procedure, which iteratively searches through all permutations of seven 



unknown blocks and 11 known blocks, and performs the sort of analysis described in Section 4.4. As 



described in Section 3.1, one may apply permutation matrices, 11 and ^P, to A and B to get 11 ^Ai> 
and ^^^Bli, and obtain at most 36 other equivalent configurations. This property was exploited to 
reduce the computations needed from 31,824 cases to about 1,500 cases. Only one matrix inverse 
completion problem was analyzed from each equivalence class. 

First, we will describe the procedures followed for a particular configuration. Then, we will list 
the pseudo-code which performed the necessary analysis for the entire problem. Our proof will be 
completed with a discussion of the results of our Mathematica procedure. 

4.5.1 A particular configuration 

We created a two-iteration, partial Grobner basis from the polynomial matrix equations resulting 
from AB and BA, along with the invertibility relations of the knowns. 
The order we used to create the Grobner basis was the following 

ki < k2 < < k4 < < kg < < ks < kg < kio < fcn 

(60) ^ Ml <C M2 ^ W3 *C M4 ^ M5 ^ Me < 

where the kj represents the j^^ known block and represents the i^^ unknown block. Inverses 
have been suppressed in our lists of knowns for clarity. Any listing of known variables should be 
accompanied by their inverses. These inverses are placed directly above, and in the same group as, 
the original variable. So, our order truly begins ki < ki^ < k2 < k2^ < . . . . 
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The output of the Grobner basis algorithm, in virtuaUy ah cases, was of the weakly essentially 



decoupled form described in Section 2.2, equations (|13|-^5|). 

To establish weakly essentially decoupled and compatibility 3-nondegeneracy, we used the output 
of the Grobner basis algorithm, which consisted solely of known indcterminates, equations (nS-na), 



or G n 5 in the language of Section 4.2. We ran the Grobner basis algorithm for one more iteration 
on these known relations, thereby creating a three iteration partial Grobner basis. We used this 
partial Grobner basis to attempt to reduce the relations which contain the unknown indcterminates. 

After applying the Grobner rules, associated with this partial Grobner basis, to the set of relations 
containing unknown indcterminates, our set of relations still had the form given in equations jl^psl). 
That is, we verified compatibility 3-nondegeneracy as given by condition (^8|). This verification was 
done by computer. This shows that the problem associated with this particular configuration is 
weakly essentially decoupled and compatibility 3-nondegenerate. 



Configuration (32) and permutations of this configuration were weakly formally backsolvable and 



compatibility 3-nondegenerate. Theorem]^ follows. 
4.5.2 Pseudo-code 

Here we give some pseudo-code with a Mathematica slant, which performs the sort of analysis 
described in the above section for all seven unknown and 11 known configurations. An essential 
part of the algorithm is the function NCMakeGroebnerBasis [polys ,k] , which creates a k iteration 
partial Grobner basis from polys. 

First, we create the relations which are implied by the invertibility of the knowns. 
inverses = NCMakeRelations [{Inv, kl,k2,k3,k4,k5,k6,k7,k8,k9,kl0,kll}] 

Next, we set the monomial order for the Grobner basis computation. This order is given in (|60|). 
SetMonomialOrder [ kl<k2<k3<k4<k5<k6<k7<k8<k9<kl0<kll«ul«u2«u3«u4«u5«u6«u7 ] 

We then generate all permutations of seven I's and eleven O's. 

permList = Permutations [ {0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1}] 

We examine all the configurations associated with the generated permutations with the 
following For[ ] loop. 

For [ i = 1, i++, i <= Length [ permList ], 



If a permutation (in the sense of Section 3.1) of this configuration was already examined, don't bother. 



If [ MemberQ[ alreadyDoneList , permList [ [i] ] ] 
Continue [] 

] 



Since no permutation of this configuration has been analyzed, we add all permutations 
of this configuration to the alreadyDoneList. 
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AppendTo [ alreadyDoneList , MakeTransf ormations [ permList [ [i] ] ] ] 



Next, convert this configuration into two Mathematica matrices. 
{A,B} = MakeSymbolicMatrices [ permList [ [i] ] ] 

Obtain the union of all relations: AB = /, BA = / and invertibility of knowns. 
relations = Union [ inverses, Flatten [ 

MatrixMultiply [ A,B ] - IdentityMatrix [3] , 

MatrixMultiply [ B,A ] - IdentityMatrix [3] ] ] 

Make a Grobner basis from the relations generated in the previous step, 
output = NCMakeGroebnerBasis [ relations, 2 ] 

Isolate the compatibility conditions on the knowns: 



G n S in the language of Section 4.2 or equations ( [13| - |15D in Section 2.2. 
polyslnKnowns = FindPolysInOnlyTheVariables [ output, 

-[kl,k2,k3,k4,k5,k6,k7,k8,k9,kl0,kll>] 



Reduce the output of the Grobner Basis Algorithm with the compatibility conditions found in the 
previous step. For our purposes, this will result in a compatibility 3-nondegenerate set of equations. 
reductionSet = NCMakeGroebnerBasis [ knownPolys, 1 ] 

output = NCReduction[ output, PolyToRule [ reductionSet ] ] 

Extract unknowns, which lie in an equation with no other unknowns, from the reduced output, 
determinedlndeterminates = PolysInOneUnknown[ output, {ul,u2,u3,u4,u5,u6,u7} ] 

Extract singleton indeterminates from the reduced output. 

s ingle Indeterminates = PolysExplicit [output , {ul ,u2 ,u3 ,u4,u5,u6 ,u7}] 

If all our unknown indeterminates lie in an equation in one unknown or are singletons, then our 

solution set is of the weakly essentially decoupled form. Otherwise, our solution set is not. 

If [Union [determinedlndeterminates , singlelndeterminates] =={ul ,u2 ,u3,u4,u5 ,u6 ,u7} , 

Print ["SUCCESSFUL"] 

] 

Else[ 

Print ["UNSUCCESSFUL"] 

] 



] (* End of For [] loop *) 
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4.5.3 End game 

The problems which were strongly undetermined did not have the formally backsolvable form. Of 
the problems which were not strongly undetermined, there were seven cases in which the output 
of the two iteration partial Grobner basis did not have the form of equations jT^ps]). For these 
seven cases, we performed the same analysis, but created a three iteration partial Grobner basis 
instead of halting the algorithm after two iterations, as was done originally. In six of these cases, the 
three iteration partial bases had the form of equations (p^p5[) and were shown to be compatibility 
3-nondegenerate. In the case associated with configuration (|3^), the 3- iteration partial Grobner 
basis did not have the essentially decoupled form, and a 4-iteration partial Grobner basis proved too 
difBcult to compute. Therefore, the result stated in the theorem follows. 

The Mathematica code, associated with the pseudo-code given above, ran for approximately 
3 days, on a Sun Ultra II with two 166Mhz processors and 1Gb of RAM. The computer was a 
departmental machine and the processes associated with these computations were therefore given 
only a portion of the total computational resources available. The same computations on a similar 
machine dedicated to this problem might take half the time. □ 

4.6 Proof of Theorem |2| 

We shall need the following lemma for our proof: 

Lemma 1 (Schur) If xi^i, xi^2, X2,i, X2.2 cire invertible block matrices of the same size, then 

Xl,l Xi,2 

2:2,1 2:2,2 

is invertible if and only if — 2:2, i2:^![ 0:1^2 + X2,2 'is invertible. 
Proof: 

/ / ^ \ ( ^I'l ^ \ f ^ ^r,l2:i,2 \ ^ [ 2:1,1 21^2 

V 2:2,125^1 ^ / V -22,l2:j'j2:i,2 + 2:2.2 J \ I J \ 2:2,1 22,2 

□ 

Proof: [Of Theorem] 

Creating a three iteration partial Grobner basis with the relations 

(61) AB = I, BA = I, and the invcrtibility of the knowns, 
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using the NCGB command NCProces£|^, yields a set of polynomials, which includes relations 

(62) z — zez + zdg + jhz — kci + jag 

(63) a-^h-'^ -a^^bjh-^ +a-^bzr^ = -d-^r^ - d-^ejh-^ + d-^ezr^ 

(64) f-^a-^ - f-^gda-^ +k-hda-^ = -k'^b-^ - f-^geb-^ + k-^eb-^ 



and relations (|42|-[47D. See Appendix 2 pages 8-10. The order used is given on page |30|, order (p5D. 
Since polynomials created through the Grobner basis algorithm are in the polynomial ideal generated 
by the original relations, the validity of relations (|6^p4[) and ( p2[^7| ) is a consequence of relations 



For us to write the relations ( p2[ - |64| ) in the form ( p5[ff0| ), we require the invertibility of {da^^ — 
eb~^) and (a^^b — d^^e). These invertibility relations are provided by the Schur lemma given above 
since the outer matrix (^^, consisting of a,b,d, and e (all knowns), is assumed to be invertible. With 
this, we can solve for z explicitly in equations ( |63[]6^ ) and write the relations (^0[pl|) defining z. 
Furthermore, we may use these definitions of z to write the relations (p2[p4) as (p£ 



The converse is again approached using a Grobner basis method. As above, the Schur complement 
formulas give the invertibility of {da~^ — eb^^) and {a^^b — d~^e), which shows that relations 
(|6^64) follow from ( p5[ - |4^ ). The question then becomes whether or not relations (^ij) arc in the 



ideal generated by polynomials ( p2[ - ^ ) and ( |6^ - |6^ ). We create a seven iteration partial Grobner 



basis 07 from the polynomials (^|7D and (|6^]64D with the NCGB command NCMakeGB, under 



the graded (length) lexicographic monomial order. One can verify that the original equations (|61 
reduce to with respect to Qt. This shows that the relations AB — I and BA = I are elements 
of the noncommutative polynomial ideal generated by the relations (|4^^7|) and (|6^|6^ and the 
invertibility of the knowns. The result follows. □ 



4.7 Discovering Theorem 121 and its proof 



In this section, we describe the process used to discover our particular theorem. Theorem 0. This 



process follows the formal notion of a strategy, rigorously developed in [HS99| 



'Appendix 2 page 7 contains the input to the NCProcess command, the "unraveled" equations (|6l|). 
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4.7.1 Addressing our problem 

In light of our goal, creating polynomials in few unknowns, we used this monomial order^ 
(65) a<5<c<d<e</<g<ft.<i<j<fc<z<u<v<w<a:<y 

and ran the Grobner basis algorithm with an iteration limit of three. 

The output of this Grobner basis computation included relations (|4^^7|), as well as ( |62[ - |6^ ). (See 
Appendix 2 pages 8-10 for the entire output of the GBA.) Thus, these relations are a consequence of 
the original relations. The necessity part of the proof is complete, modulo a bit of Schur complement 



beautification done in Section 4.7.4 



4.7.2 Converse: a smaller basis 

It is true that the original relations ( |6l| ) are members of the ideal generated by the long and ugly 
relations taking up pages 8-10 of Appendix 2 (The partiality of a Grobner basis at some iteration 
is only in respect to its reduction properties, and not the ideal generated by these relations.). We 
could have written these down instead of equations (|3^p7[), our final conclusion, and stopped, but 
we would prefer to have a more concise set of relations which imply the original relations. In other 
words, we would like to have a smaller basis for this ideal. 

The computer commands in NCGB have the ability to simplify the basis in the manner above, 
in the same step as generating it, by setting certain options. However, we did not have the com- 
puting power, or perhaps the patience, to isolate the few relations on z given above ( |62[]6^ ) using 
this method, under the original order. To this end, the monomial order was changed to graded 
lexicographic. In NCGB notation, we replaced all of the <C's with <'s. The graded lexicographic 
order computations are often of much less computational complexity, since monomials usually must 
be merely checked for number of elements. When our original order was imposed on the small 
basis algorithm, the two iteration application did not complete after several days running on a Sun 
SPARCstation-4 computer, while under the graded lexicographic order the algorithm finished in a 
few minutes. 

We tried several different sequences, of which most gave unsatisfactory results. The bases found 
were not small enough in these cases. An acceptable small basis obtained through this procedure 
consisted of the invertibility relations on the knowns, the relations which give the unknowns other 
than z in terms of z ( ^2[^7|) , and relations concerning z and the knowns (|6^|64|). The computer 
work associated with this is given in Appendix 3. 



* Inverses have been suppressed in our lists of knowns for clarity. Any listing of known variables should be 
accompanied by their inverses. These inverses are placed directly above and in the same group as the original 
variable. So our order truly begins a < < b < < . . . . 
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4.7.3 Confirmation 



To confirm that these relations (62-BJ, [4^ - |4^ , and invertibihty of the knowns) imply the original 
relations, we created a noncommutative partial Grdbncr basis from these relations and reduced the 
original relations with this partial Grobner basis. The original relations all reduced to 0. Thus, it was 
shown that the original relations were elements of the ideal generated by the relations given above. 
Interestingly enough, a five iteration partial Grobner basis did not reduce the original relations 
although a seven iteration partial Grobner basis did. (See Appendix 4.) The order used for this 
computation was again the graded lexicographic. 

4.7.4 Beautification with Sciiur Complements 



Equations (|63[|6j) are especially appreciated, because they are linear in one unknown variable z. 
A more satisfying situation, though, would be to have an expression for z entirely in terms of the 
knowns, a singleton equation. This may be accomplished by assuming the invertibihty of 

{prH-d-^e) and {doT^ - eh-^) 

in equations (|6^Q). By the Schur Lemma |l|, this is equivalent to the invertibihty of the outer 
matrix ( ^ ^ I i since all entries of this matrix are themselves invertible. At the outset of our 



d e 

investigations, we had no reason to assume this more restrictive condition. It was only after realizing 
the utility of this assumption that we added it to our conditions. 



With this, z is given explicitly by the equations (41 -401) and each of these must satisfy the 
quadratic ([6^). Hence, the equations (35-3^) on the knowns are a necessary and sufficient set of 
conditions for AB — I and BA — I. 

5 Conclusion 

In this article, we have investigated the use of noncommutative symbolic algebra software in the 
analysis of partially prescribed inverse matrix completion problems. We described a method for 
solving such problems with a computer. We have shown that the solutions to all 3x3 block inverse 
matrix completion problems with seven unknown and 11 known blocks are of a relatively simple form. 
We presented one particular theorem, and showed how it can be massaged into a more palatable 
form by making some mild assumptions on the prescribed (known) blocks. 

Finally, the author would like to express appreciation for the effort put forth by the anonymous 
referee. His or her careful reading, criticism, and editing have greatly improved this paper. 
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Appendices: Using noncommutative Grobner bases in solving 
partially prescribed matrix inverse completion problems 



February 1, 2008 



We intend these appendices to appear on the internet for those 
interested and not be included in the pubhcation. 







Abstract 

We investigate the use of noncommutative Grobner bases in solving partially prescribed 
matrix inverse completion problems. The types of problems considered here are similar to those 
in [B JEW] . There the authors gave necessary and sufficient conditions for the solution of a two 
by two block matrix completion problem. Our approach is quite different from theirs and relies 
on symbolic computer algebra. 

Here we describe a general method by which all block matrix completion problems of this 
type may be analyzed if sufficient computational power is available. We also demonstrate our 
method with an analysis of all three by three block matrix inverse completion problems with 
eleven known blocks and seven unknown. We discover that the solutions to all such problems 
are of a relatively simple form. 

We then perform a more detailed analysis of a particular problem from the 31,824 three 
by three block matrix completion problems with eleven known blocks and seven unknown. A 
solution to this problem of the form derived in [BJLW] is presented. 

Not only do we give a proof of our detailed result, but we describe the strategy used in 
discovering our theorem and proof, since it is somewhat unusual for these types of problems. 
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1 Appendix 1 - Mathematica implementation of pseudo code 
given for Theorem 1 



(* First we set up needed variables and lists *) 

onePerms = Permutations [{0,0, 0,0,0, 0,0, 0,0,0,0, 1 , 1 , 1 , 1 , 1 , 1 , 1}] ; 

vars = { kl,k2,k3,k4,k5,k6,k7,k8,k9,kl0,kll,ul,u2,u3,u4,u5,u6,u7 }; 

SetNonCommutative [ vars ]; 

SetNonCommutative [ Inv [kl] , Inv [k2] , Inv [k3] , Inv [k4] , Inv [k5] , Inv [k6] , Inv [k7] , 
Inv [k8] , Inv [k9] , Inv [klO] , Inv [kll] ] ; 

inversesKeep = 

{-l+kl**Inv[kl] ,-l+klO**Inv[klO] ,-l+kll**Inv [kll] ,-l+k2**Inv[k2] , 
-l+k3**Inv[k3] , -l+k4**Inv [k4] ,-l+k5**Inv[k5] , -H-k6**Inv [k6] , 
-l+k7**lnv[k7] , -l+k8**lnv [k8] , -l+k9**lnv [k9] , -1+lnv [kl] **kl , 
-1+Inv [klO] **klO , -1+Inv [kll] **kll , -1+Inv [k2] **k2 , -1+Inv [k3] **k3 , 
-1+Inv [k4] **k4 , -1+Inv [k5] **k5 , -1+Inv [k6] **k6 , -1+Inv [k7] **k7 , 
-1+Inv [k8]**k8, -1+Inv [k9]**k9 >; 

permMtcs = Permutations [ IdentityMatrix [3] ]; 

invPermMtcs = Map[ Inverse, permMtcs ]; 

knownVars = {kl ,k2,k3,k4,k5,k6,k7,k8,k9,kl0,kll}; 
unkVars = {ul,u2,u3,u4,u5,u6,u7 }; 

alreadyDone = {}■ ; 

(* Here we give some necessary functions *) 

(* Make symbolic matrix from I's notation *) 

MakeSymbMatrix [ onesList_List ] := Module [{idx,unkldx=l ,knldx=l ,newMtx={} }, 

For [idx=l , idx<=Length [onesList] , idx++ , 

If [ onesList [[idx] ] == 1, 

AppendTo [ newMtx , unkVar s [ [unkldx] ] ] ; 
unkldx++ ; , 

AppendTo [ newMtx , knownVars [ [knidx] ] ] ; 
knldx++; 

]; 

]; 

Return [ newMtx ] ; 

]; 

(* Ask if a permutation of a matrix is in the alreadyDone list *) 
PermMemberQ [ onesList_List ] := 
Module [ {permList={}, A,B>, 

A = Partition [onesList, 9 ] [[1]]; 
B = Partition [onesList, 9 ][[2]]; 
A = Partition [A, 3] ; 
B = Partition [B, 3] ; 

Apermed = Map[ A.# & , permMtcs] ; 

Aperm2 = Flatten[Table [Map [#. Apermed [[permEntry]]&, permMtcs], 



2 



{permEntry , 1 , Length [Apermed] }],!]; 

Bpermed = Map[ #.B & , invPermMtcs ]; 

Bperin2 = Flatten [ Table [Map [ Bpermed [ [permEntry] ]. # &, invPermMtcs ], 
{ permEntry , 1 , Length [Bpermed] }],!]; 

For [index =1 , index<=Length [Aperm2] , index++ , 

AppendTo [permList , Flatten [-[Aperm2 [ [index] ] , Bperm2 [ [index] ] }] ] ; 

]; 

If [ Intersection[ alreadyDone, permList ] === {}, 

Return [ False ] ; , 
Return [ True ] ; 

]; 

]; 

(* Find the unknown indeterminates which lie in equations in one unknown *) 

FindDetermined [ aList_List , currentUnknowns_List ,opts Rule] := 

Module [{i , len, item, VEirs ,kn,unk,n,alllengths , rules , determList , 
allvars , allind , j , total vars , outputToFile , f ileMame} , 

Clear [relations] ; 
determList = {}; 

rules = Union [ExpandNonCommutativeMultiply [aList] ] ; 
rules = PolyToRule [rules] ; 
rules = Union [rules] ; 
len = Length [rules] ; 

Do [item = rules [[i]]; 

If [Mot [item===0] , 

vars = Grablndeterminants [item] ; 

If [ Length [ Union [vars] ] == 1 , 

AppendTo [determList , vars [ [1] ] ] ; 

]; 

]; 

,{i,l,len}]; (* End Do[] loop *) 
Return [ determList ] ; 

]; 

(* Find the indeterminates which are singletons *) 

FindSingletons [ aList _List , currentUnknowns_List ,opts Rule] := 

Module [{i, len, item, vars, kn,unk,n,alllengths, rules, singletonList , 
allvars , allind , j , total vars , outputToFile , f ileMame} , 

Clear [relations] ; 
singletonList = {}; 

rules = Union [ExpandNonCommutativeMultiply [aList] ] ; 
rules = PolyToRule [rules] ; 
rules = Union [rules] ; 
len = Length [rules] ; 

Do [item = rules [[i]]; 

If [Mot [item===0] , 

vars = Grablndeterminants [item] ; 
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(* If the head is not NCM it's a singleton !! *) 
If [Head [item [[!]]]=! =NonCommutativeMultiply , 
AppendTo [singletonList , iteiii[[l]] ]; 

]; 

]; 

,{i,l,len>]; (* End Do [] loop *) 

Return [Intersection [ singletonList, currentUnknowns ] 

]; 

]; 

<<Extra . TeXForm ; 

WriteMatrixTeX [ A_List, B_List, currNum_ ] := Module [{}, 

OpenWriteForTeX["thms/mtcs"<>ToString[currNum] ] ; 
ExpressionToTeXFile["thms/intcs"<>ToString[currNum]<>" .tex" , 
OutputAMatrix [A] ] ; 

ExpressionToTeXFile["thms/intcs"<>ToString[currNum]<>" .tex" , 
OutputAMatrix [B] ] ; 

WriteString["thms/mtcs"<>ToString[currMum]<>" .tex", 
"\\end{document}-" ] ; 

Close["thms/iiitcs"<>ToString[currNuin]<>" .tex"] ; 

]; 

OpenWrite ["AnswersForCompleteMtx"] ; 
undetermList = 

{683 , 684 ,689,695,706,708,710,719, 725 , 729 , 748 , 749 , 752 , 755 ,760,769, 779 , 784 , 1656 , 
1675 , 1701 , 2580 , 2599 , 2625 , 4296 , 4315 , 4341 , 5957 , 5958 , 5963 , 5969 , 5980 , 5982 , 5984 , 
5993,5999,6003,6022,6023,6026,6029,6034,6043,6053,6058,6090,6109,6135,6174, 
6193 , 6219 , 6294 , 6313 , 6339 , 7299 , 7318 , 7344 , 7486 , 7487 , 7491 , 7498 , 7499 , 7500 , 7501 , 
7520 , 7521 , 7540 , 7556 , 7575 , 7606 , 7607 , 7608 , 7609 , 7610, 7611 , 7738 , 7757 , 7821 , 7948 , 
7967 , 8031 , 9070 , 9089 , 9153 , 10357 , 10376 , 10440 , 11244 , 11263 , 1 1289 , 12304 , 12323 , 
12349 , 13409 , 13410 , 13414 , 13421 , 13422 , 13423 , 13424 , 13443 , 13444 , 13463 , 13479 , 
13498 , 13529 , 13530 , 13531 , 13532 , 13533 , 13534 , 13535 , 13554 , 13618 , 13829 , 13848 , 
13912 , 14159 , 14178 , 14242 , 15941 , 15960 , 16024, 16249 , 16268 , 16294 , 17364, 17383 , 
17447,17943,17962,18026,20312,20331,20357,22914,22915,22919,22926,22927, 
22928 , 22929 , 22948 , 22949 , 22968 , 22984 , 23003 , 23034 , 23035 , 23036 , 23037 , 23038 , 
23039 , 23040 , 23059 , 23123 , 23124 , 23143 , 23207 , 23574 , 23593 , 23657 , 24069 , 24088 , 
24152 , 24257 , 24276 , 24302 , 25372 , 25391 , 25455 , 26786 , 26805 , 26869 , 28954 , 28973 , 
29037 , 29789 , 29808 , 29872 , 30476 , 30477 , 30482 , 30488 , 30499 , 30501 , 30503 , 30512 , 
30518 , 30522 , 30541 , 30542 , 30545 , 30548 , 30553 , 30562 , 30572 , 30577 , 30609 , 30628 , 
30654 , 30693 , 30712 , 30738 , 30813 , 30832 , 30858 , 30978 , 30997 , 31023 , 31 198 , 31217 , 
31243, 31484, 31503, 31529> ; 



(* Here is the MAIN LOOP *) 

For [ curr Index = 1 , currlndex <= Length [onePerms] , 
curr Index ++, 

Write ["AnswersForCompleteMtx" , "permutation*" , 
currlndex ] ; 

If [ PermMemberQ [ onePerms [ [currlndex] ] ] , 
Print [ "Found One Already Done MM " ] ; 
Write ["AnswersForCompleteMtx" , 
" Found one already done", 
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onePerms [ [curr Index] ] ] ; 
Continue [] ; 

]; 

(* Only add to alreadyDone list if a perm is not in it *) 
AppendToE alreadyDone, 

onePerms [ [currlndex] ] ] ; 

If [ MemberQ [undetermList , currlndex ] , 
Write ["AnswersForCompleteMtx" , 

"This problem is of UNDETERMINED FORM. 

onePerms [ [currlndex] ] ] ; 
Continue [] ; 

]; 

alreadyDone = Union [ alreadyDone ] ; 
newMtcs = MakeSymbMatrix [onePerms [ [currlndex] ] ]; 
mtrcs = Partition [ newMtcs ,9 ] ; 

matrixA = Partition [ mtrcs [[1]] ,3 ] ; 
matrixB = Partition [ mtrcs [[2]] ,3 ] ; 

WriteMatrixTeX[matrixA, matrixB, currlndex ]; 

oneway = MatMult [matrixA, matrixB] - IdentityMatrix [3] ; 
otherway = MatMult [matrixB, matrixA] - IdentityMatrix [3] ; 
start = Flatten [{ oneway, otherway }] ; 
inverses = inversesKeep; 
start=Join [start, inverses] ; 

ClearMonomialOrderAll [] ; 

SetMonomialOrder [ {{kl , Inv [kl] , k2 , Inv [k2] , k3 , Inv [k3] , k4 , Inv [k4] , 
k5 , Inv [k5] , k6 , Inv [k6] , k7 , Inv [k7] , k8 , Inv [k8] , k9 , Inv [k9] , 
klO , Inv [klO] ,kll , Inv [kll] } , {ul} , {u2> , {u3} , {u4} , {u5} , {u6} , {u7}}] ; 

f ileName = "thms/WoerdOutput-"<>ToString [currlndex] ; 

(* This function creates a Noncommutative Groebner basis and 
a TeX file describing the output *) 
answer = NCProcess [start ,2,2, 1 , 1 , 

f ileName ,RR->True, SB->True, SBByCat->False, NCCV->False 

] ; 

If [ Complement [unkVars , Union [ 

FindDetermined [answer [[3]] , unkVars] , 
FindSingletons [ answer [[3]], unkVars] ] ] =!= {} , 

Write ["AnswersForCompleteMtx" , "Didn't WORK correctly.", 

onePerms [ [currlndex] ] ] ; , 
(* else *) 

Write ["AnswersForCompleteMtx" , "Worked fine.", 
onePerms [ [currlndex] ] ] ; 

] ; (* end If [] *) 
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Open ["thms/outRels"<>ToString [curr Index] <>" .m"] ; 
Write ["tlims/outRels"<>ToString [curr Index] <>" .m" , 

answer [ [3] ] ] ; 
Close ["thms/outRels"<>ToString[currIndex] <>" .m"] ; 

] ; (* end for *) 

Close [ "AnswersForCompleteMtx"] ; 



(^^^!ic^^^^^^ Redefine SmallBasisC ] for decoupled analysis *****) 
SmallBasis [input_List,keepListInverses_List , iterationCount_Integer] := 
Module [ O, 

singlePolys = FindSinglePolys [ result, -[ ul,u2,u3,u4, 
u5,u6,u7 > ] ; 

result = Complement [ result , singlePolys ] ; 
keepList = inversesKeep ; 
ClearMonomialOrderAll [] ; 

SetMonomialOrderC { kl,Inv[kl] ,k2,Inv[k2] ,k3,Inv[k3] ,k4, 
Inv [k4] , k5 , Inv [k5] , k6 , Inv [k6] , k7 , Inv [k7] , 
k8,Inv[k8] ,k9,Inv[k9] ,klO,Inv[klO] ,kll , Inv [kll] , 
ul,u2,u3,u4,u5,u6,u7}] ; 

knownPolys = FindKnownPolys [ result, {. ul,u2,u3,u4, 
u5,u6,u7 } ] ; 

(* Here we order the polys since order matters for the Small Basis Alg *) 
result = Flatten [ { knownPolys, 

FindDeterminedPolys [ result, { ul,u2,u3,u4, 
u5,u6,u7 } ] , 

Complement [ result. Union [ knownPolys , 

FindDeterminedPolys [ result, { ul,u2,u3,u4, 
u5,u6,u7 } ] 
] 

] > 

]; 

keepList = Union [ keepList, knownPolys ]; 

wRules = PolyToRule[ NCMakeGB[ keepList, 1 ] ]; 

result = Reduction [ Complement [ result , keepList ], wRules ]; 

result = Union [ result, keepList, singlePolys ]; 

(* Restore normalcy for Regular Output to look nice *) 
result = Union [ result, singlePolys, keepList ]; 

ClearMonomialOrderAll [] ; 

SetMonomialOrder [ {kl,Inv[kl] ,k2,Inv[k2] ,k3,Inv[k3] ,k4, 
Inv [k4] , k5 , Inv [k5] , 

k6,Inv[k6] ,k7,Inv[k7] ,k8,Inv[k8] ,k9,Inv[k9] ,klO,Inv[klO] , 
kll , Inv [kll] } , {ul} , {u2} , {u3} , {u4} , {u5} , {u6} , {u7} ] ; 
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]; 



2 Appendix 2 - The First Run for Theorem 2 



Input = 

— I + ax + bj + th 
a f + bk + ty 

ag + bz + ti 
ch + ux + V j 

— 1 + cy + uf + vk 
ci + u g + V z 

dx + ej + wh 
d f + ek + wy 

— 1+dg + ez + wi 

— 1 + fu + gd + xa 
f c + g w + xt 

f V + g e + xb 
ha + i d + y u 

— 1 + ht + iw + yc 
hb + ie + yv 

j a + ku + zd 
j t + k c + zw 

— 1+jb + kv + ze 
aa~^ == 1 

a^^ a == 1 
66-1 I 
6-1 b == 1 
cc-i == 1 
c-^c== 1 
dd-^ == 1 
d-^d== 1 
ee-i == 1 
e == 1 
==1 

gg-^ == 1 

/i/i-i == 1 
h-^h== 1 
ii-i == 1 
i-M == 1 

J == 1 

== 1 

kr^ k == 1 

File Name = MatrixInvorscAnswcr-3 
NCMakcGB Iterations = 2 
NCMakeGB on Digested Iterations = 3 
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SmallBasis Iterations = 3 
SmallBasis on Knowns Iterations = 4 
Deselect—* {} 

UserSelect^ {} 

RR^Tnic 

RRByCat^False 

SB^False 

SBByCat^False 

DegrccCap^l2 

DegrccSurnCap^SO 

DegreeCapSB^ia 

DegreeSumCapSB^Sl 

NCCV^Falsc 

THE ORDER IS NOW THE FOLLOWING: 

a < < 6 < < c < < d < < e < < / < < g < g~'^ < h < h~'^ < z < < 

j < j^^ < k < 1:^^- -C ~ <r // < .r < f < II < r < ir 



YOUR SESSION HAS DIGESTED 
^^^^^^^^^^^^^ THE FOLLOWING RELATIONS ^^^^^^^^^^^^ 

THE FOLLOWING VARIABLES HAVE BEEN SOLVED FOR: 
{t, u, V, w, X, y} 

The corresponding rules are the following: 
t — > —lagi~^ — bzi~^ 

—lk~^ja — k~^zd 

V — > — jb — k~^ z e 

w — dgi~^ — ezi~^ 

X + f k~^ j — gda~^ + f k~^ zda~^ 

y ^ k~^ j a f + k~^ j bk + k~^ zd f + k~^ zek 

The expressions with unknown variables {} 

and knowns {a, b, c, d, e, /, 5, h, i, j, k, a~^, b~^, c~^, d~^, e~^, f~^, g~^, h~^, k~^} 

aa~^ — > 1 

66-1 ^ 1 

cc-i ^ 1 

dd-^ 1 

ee-i 1 

99-'^ I 
hh-'^ 1 
ii-^ 1 

fcfc-i ^ 1 

a — > 1 
b-^b^l 
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ft-i h 1 

USER CREATIONS APPEAR BELOW 

SOME RELATIONS WHICH APPEAR BELOW 
^^^^^^^^^^^^^^M MAY BE UNDIGESTED ^^^^^^^^^^^^^^M 

THE FOLLOWING VARIABLES HAVE NOT BEEN SOLVED FOR: 

{a, h, c, d, e, /, g, h, i, j, k, z, a'^, b''^, c'^, d'^, e'^, /"S g-'^, h'^, r\ fc-^} 

The expressions with unknown variables {z} 

and knowns {a, b, c, d, e, f, g, h, i, j, k, a~^, b~^, c~^, d~^, e~^, f~^, h~^, k~^} 
zez^z — jag — jbz + kci — zdg 

d-^ezi-'^ ^ a-^h-^ -a-^bjh-^ +a-^bzi-'^ +d-^ejh-^ 

fc-i zeb-^^ /-I + fc-i - f-'^gda'^ +f~'^9e 6"^ + fc-^ zda'^ 

iezi~^c~^ —lhagi~^c~^ — hbzi~^ — idgi~^c~^ + c~^k~^jaf + c~^k~^jbk + 

k~^ zdf + k~^ zek 

k-^zda-^b-^ -Ifc-i - f-^ ge- f-^ a-^b + k-^ ze + f-^ gda-H 

k~^ zekc ^ chagi~^ + chb zi~^ + cidgi~^ + cie zi~^ — k~^ j a f c — k~^ j bkc — k~^ zd f c 

afk~^zda~^ —Ibj — afk~^j + agda~^ + agi~^ h + bzi~^ h 

dfk~^zda~^ —lda~^ — ej — h — df k~^ j + dgda~^ + dgi~^ h + ezi~^ h 

fk-'^zda-'^h-^ gi-^ - f k-'^ j h-'^ + g da-^ h-^ - a-^ b j h-'^ + a-^ b zi-'^ 

i e z h a ^ h a + i d — h a g d + h b j a — i d g d + i e j a + h a f k~^ j a + h a f k~^ z d — h a g h a — 
hb zi~^ ha + i d f k~^ j a + i d f k~^ z d — i d gi~^ ha 

ie zi~^ hb ^ —lha f k~^ — hag e + hbj b — id f k~^ — idg e + iej b + ha f k~^ j b + ha f k~^ ze — 
hagi~^ hb — hbzi~^ hb + id f k~^ j b + id f k~^ ze — idgi~^ hb 

kciezi~^ ^ j af c + jbkc + zdf c + zekc — kchagi~^ — kchbzi~^ — kcidgi~^ 

z d f k~^ z d ^ —Izd+jagd—j bj a+kcha+z dg d—z e j a—j a f kr^ j a—j a f k~^ z d—z d f k~^ j a 

z d f k~^ z e j b + j a f k^^ + jage — jbjb + kchb + zdfk^^ + zdg e — z e j b — j a f k~^ j b — 
j af k~^ ze — zdfk~^jb 

k~^ z d f k~^ z — 1 fc"i z + chadr^ + kr^ j ag + k^^ zdg — k~^ j a f k~^ z — k~^ jbj ad~^ — 
k~^ zej ad~^ — k^^ j a f k~^ j ad~^ — k~^ z d f k~^ j a d~^ 

bzi~^ hb zi~^ ^ a f c — b zi~^ — agdgi~^ — agezi~^ + bj agi~^ + bj b zi~^ — agi~^ hagi~^ — 
agi~^ hbzi~^ — bzi~^hagi~^ 

e zi~^ hb zi~^ dfc+ dgi~^ — dgdgi~^ — dgezi~^ + ejagi~^ + ejbzi~^ + ha g i~^ + 
i~^ hbzi~^ — dgi~^ hagi~^ — dgi~^ hbzi~^ — ezi~^ hagi~^ 
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zi ^ hbzi ^ c ^ — > —Izi ^ c ^ + b ^ a f + j agi ^ c ^ + jbzi ^ c ^ — zi ^ hagi ^ c ^ — 
b~^ agdgi~^ — b~^ agezi~^ —b~^agi~^ hagi~^ — b~^ agi~^ hbzi~^ 

z d f j be~^ ^ z — chad~^ + chbe~^ + j be^^ + k^^ j a f k^^ e^^ + kr^ j b j a dr^ — 
k~^jbjbe~^ + k~^zdfk~^e~^ + k~^zejad~^ — k~^ zejbe~^ + k~^ j a f k~^ j ad~^ — 
j a f j be~^ + k~^ zdfk~^jad~^ 

dgi~^ hb zi~^ zi~^c~^ — b~^af + e~^df + er^dgi~^c~^ + b~^ a g d g + 

b~^ a g e zi~^ —e^^ d g d g —e"^ d g e z c~^+e"^ i^^ ha gi~^ c~^+e~^ hbzi~^ c~^ + 

b^^ a g i^^ h a g i^^ c^^ + b^^ a g i^^ hb zi^^ c^^ — e^^ d g h a g i^^ c^^ 

hbzi~^ k~^ z d ^ d + ha — dgi~^ k~^ j a — dgi~^ k~^ z d — e zi~^ k~^ j a — 
ezi~^ k~^ z d — hagi~^ k~^ j a — hagi~^ k~^ z d — hbzi~^ k~^ j a 

hb zi~^ k~^ z e — *■ e + hb + dgi~^c~^k~^ + ezi~^c~^k~^ — 
dgi~^ k~^ j b — dgi~^ k~^ z e — e zi~^ j b — e zi~^ k~^ z e + i~^ hagir^ c^^ k~^ + 

hb zi~^ k~^ — i^^ ha g k~^ j b — ha g i^^ c^^ k^^ z e — hb zi^^ c^^ k^^ j b 

k~^ zdf k~^ j ad~^ e —Ichb — k~^ j b — k~^ze + chad~^e — k~^jafk~^ + k~^jbjb — 
k~^ zd f k~^ + k~^ zej b + k~^ j a f k~^ j b — k~^ j bj ad~^ e + k~^ zd f k~^ j b — k~^ zej ad~^ e — 
k~^ j a f k~^ j ad~^ e 

k~^ z e j a f k~^ z d —Ichagd + chbja + k~^jbja + k~^zdgd + chafk~^ja + 

c h a f k~^ z d — k~^ j a f c h a + k~^ j a f k~^ j a + k~^ j a f k~^ z d — k~^ j a g d g d + k~^ j a g e j a + 
k~^ j bj a g d — kr^jbjbja — kr^ z d f cha — kr^zdgdgd + k~^ z d g e j a + k~^ z e j a g d — 
k~^ z ej b j a + k~^ j agd f k^^ j a + k~^ j ag d f k^^ z d — k~^ j bj a f k~^ j a — k~^ j bj a f k~^ z d + 
k~^ z d g d f k^^ j a + k~^ z d g d f k~^ z d — k~^ z e j a f k~^ j a 

k~^ z ej a f k~^ z e —Ichb — k~^jb — chafk~^ — chage + chbjb — k~^jafk~^ — 
k~^ j ag e+2k~^ j bj b+k~^ z ej b+cha f k~^ j b+cha f k~^ ze — k~^ j a f chb+k~^ j a f k~^ j b+ 
k~^ j a f k~^ ze — k~^ j agdf k~^ 

— k~^jagdge + k~^jagejb + k~^ jbj af k~^ + k~^ jbjage — k~^ jbjbjb — k~^ zdf chb — 

k~^zdgdfk~^ — k~^zdgdge + k~^zdgejb + k~^zejafk~^ + k~^zejage — 
k~^ z e j b j b + k~^ j a g d f k~^ j b + k~^ j a g d f k~^ z e — k^^ j b j a f k~^ j b — k^^ j b j a f k^^ z e + 
k~^ zdg d f k~^ j b + k~^ zdg d f k~^ ze — k~^ z ej a f k~^ j b 

k~^ z ej bj b zi~^ ^ cha f c + chagi~^ — k~^ z d f c — k~^ zdgi~^ — cha g dg — chage zi~^ + 
chbjagi"^ + chb jb zi~^ + k^^jagdfc + 

k~^ j ag dgi^^ ~k^^ j bj a f c+k^^ j bj b zi~^ + k~^ z dg d f c+2 k~^ z d g d g + k~^ zdgezi~^ — 
k~^ z e j a f c — k^^zejagi^^ — k^^ j a f chagi~^ — k^^ j a f chb zi"^ — k~^ j agdgdgi~^ — 
k~^ j a g dg e zi~^ + kr^ j a g e j a gi^^ + kr^ j a g e j b zi~^ + kr^ j bj a g d gi^^ +kr^ j bj a g e zi~^ — 
k~^ jbjbjagi~^ — kr^ j bjbjbzi'^ — kr^zdfchagi~^ — k^^zdfchbzi^^ — kr^zdgdgdgi~^ — 
k~^ zdgdgezi~^+k~^ z dg ej agi~^+k~^ z dg ej b zi~^+k~^ z ej ag dgi~^+k~^ zejagezi~^ — 
k^^ zej bjagi~^ 

The time for preliminaries was 0:00:01 

The time for NCMakeGB 1 was 0:00:00 

The time for Remove Redundant 1 was 0:00:00 

The time for the main NCMakeGB was 0:00:05 

The time for Remove Redundant 2 was 0:00:00 

The time for reducing unknowns was 0:00:01 

The time for clean up basis was 0:00:02 

The time for SmallBasis was 0:00:01 

The time for CreateCategories was 0:01:07 

The time for NCCV was 0:00:00 

The time for RegularOutput was 0:00:38 

The time for everything so far was 0:01:58 
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3 Appendix 3 - Find A Smaller Basis 



In[l] := Input: 

SetNonCommutative [a,b,c,d,e,f,g,h,i,w,x,y,z,j,u,t,k,v, 
Inv [ a] , Inv [b] , Inv [ c] , Inv [d] , Inv [e] , Inv [f], Inv [g] , Inv [h] , 
Inv[i] , Inv[ j] , Inv[k] ]; 

(*Here are the relations which we like and would like to retain*) 
hopepolys = {-1 + a**Inv[a], -1 + b**Inv[b], -1 + c**Inv[c], 
-1 + d**Inv[d], -1 + e**Inv[e], -1 + f**Inv[f], 

-1 + g**Inv[g], -1 + h**Inv[h], -1 + i**Inv[i], 

-1 + j**Inv[j], 3-1 + k**Inv[k], -1 + Inv[a]**a, -1 + Inv[b]**b, -1 + Inv[c]**c, -1 + Inv[d]**d, -1 + Inv[e]**e, 
-1 + Inv[f]**f, -1 + Inv[g]**g, -1 + Inv[h]**h, -1 + Inv[i]**i, 
-1 + Inv[j]**j, -1 + Inv[k]**k, t + a**g**Inv[i] + b**z**Inv[i], 
w - Inv[i] + d**g**Inv[i] + c**z**Inv[i], 
-z + j**a**g + j**b**z - k**c**i + z**d**g + z**c**z, 
u + Inv[k]**j**a + Inv[k]**z**d, v - Inv[k] + Inv[k]**j**b + Inv[k]**z**e, 
X - Inv[a] - f**Inv[k]**j + g**d**Inv[a] - f**Inv[k]**z**d**Inv[a], 
y - Inv[c]**Inv[k]**j**a**f - Inv[c]**Inv[k]**j**b**k - 
Inv[c]**Inv[k]**z**d**f - Inv[c]**Inv[k]**z**e**k 

} 

(* Here are the relations resulting from the original "discovery" run *) 

allpolys = {-1 + a**Inv[a], -1 + b**Inv[b], -1 + c**Inv[c], -1 + d**Inv[d], 
-1 + c**Inv[c], -1 + f**Inv[f], -1 + g**Inv[g], -1 + h**Inv[h], 
-1 + i**Inv[i], -1 + j**Inv[j], -1 + k**Inv[k], -1 + Inv[a]**a, 
-1 + Inv[b]**b, -1 + Inv[c]**c, -1 + Inv[d]**d, -1 + Inv[e]**e, 
-1 + Inv[f]**f, -1 + Inv[g]**g, -1 + Inv[h]**h, -1 + Inv[i]**i, 
-1 + Inv[j]**j, -1 + Inv[k]**k, t + a**g**Inv[i] + b**z**Inv[i], 
w - Inv[i] + d**g**Inv[i] + c**z**Inv[i], 
-z + j**a**g + j**b**z - k**c**i + z**d**g + z**e**z, 
u + Inv[k]**j**a + Inv[k]**z**d, v - Inv[k] + Inv[k]**j**b + Inv[k]**z**e, 
-Inv[a]**Inv[h] - Inv[d]**Inv[i] + Inv[a]**b**j**Inv[h] - 
Inv[a]**b**z**Inv[i] - Inv[d]**c**j**Inv[li] + Inv[d]**e**z**Inv[i], 
-Inv[f]**Inv[a] - Inv[k]**Inv[b] + Inv[f]**g**d**Inv[a] - 
Inv[f]**g**e**Inv[b] - Inv[k]**z**d**Inv[a] + Inv[k]**z**e**Inv[b], 
X- Inv[a] - f**Inv[k]**j + g**d**Inv[a] - f**Inv[k]**z**d**Inv[a], 
y - Inv[c]**Inv[k]**j**a**f - Inv[c]**Inv[k]**j**b**k - 
Inv[c]**Inv[k]**z**d**f - Inv[c]**Inv[k]**z**c**k, 
Inv[k] + Inv[f]**g**e + Inv[f]**Inv[a]**b - Inv[k]**z**e - 
Inv[f]**g**d**Inv[a]**b + Inv[k]**z**d**Inv[a]**b, 
h**a**g**Inv[i]**Inv[c] + h**b**z**Inv[i]**Inv[c] + 
i**d**g**Inv[i]**Inv[c] + i**e**z**Inv[i]**Inv[c] - 
Inv[c]**Inv[k]**j**a**f - Inv[c]**Inv[k]**j**b**k - 
Inv[c]**Inv[k]**z**d**f - Inv[c]**Inv[k]**z**e**k, 
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In[2]:= ~c**h**a**g**Inv [i] - c**h**b**z**Inv [ i ] - c**i**d**g**Inv [i] - 
c**i**g**z**Inv [i] + Inv [k] ** j**a**f **c + Inv [k] ** j**b**k**c + 
Inv [k] **z**cl**f **c + Inv [k] **z**e**k**c, 

b**j + a**f **Inv [k] ** j - a**g**d**Inv [a] - a* *g**Inv [ i ] **h - 
b**z**Inv [i] **h + a**f **Inv [k] **z**d**Inv [a] , 

d**Inv[a] + e**j + Inv[i]**h + d* * f * * Inv [k] * * j - d* *g* *d* * Inv [ a ] - 
d**g**Inv [i] **h - e**z**Inv [i] **h + d**f **Inv [k] **z**d**Inv [a] , 
-g**Inv[i] + f**Inv[k] ** j**Inv[h] - g**d**Inv [a] **Inv [h] + 
Inv [a] **b** j **Inv [h] - Inv [ a ] **b**z**Inv [ i ] + 
f**Inv[k]**z**d**Inv[a]**Inv[h], 

-h**a - i**d + h**a**g**d - h**b**j**a + i**d**g**d - i**e**j**a - 
h**a**f**Inv[k] **j**a - h**a**f **Inv [k] **z**d + h**a**g**Inv [i] **h**a + 
h**b**z**Inv[i] **h**a - i**d**f **Inv [k] ** j **a - i**d**f **Inv [k] **z**d + 
i**d**g**Inv [i] **h**a + i**e**z**Inv [i] **h**a, 

h**a**f**Inv[k] + h**a**g**e - h**b**j**b + i**d**f **Inv [k] + 
i**d**g**e - i**e**j**b - h**a**f **Inv [k] ** j**b - 

h**a**f**Inv[k] **z**e + h**a**g**Inv [i] **h**b + h**b**z**Inv [i] **h**b - 
i**d**f**Inv[k] **j**b - i**d**f **Inv [k] **z**e + i**d**g**Inv [i] **h**b + 
i**e**z**Inv [i] **h**b, -j**a**f**c - j**b**k**c - z**d**f**c - 
z**e**k**c + k**c**h**a**g**Inv [i] + k**c**h**b**z**Inv [ i ] + 

]^**Q**jL**^**g** jj^-^^ [jl] + ]^ + + g + + -j^ + + g + + 2**Inv [ i ] , 

z**d - j**a**g**d + j**b**j**a - k**c**h**a - z**d**g**d + z**e**j**a + 
j**a**f**Inv[k] ** j**a + j**a**f **Inv [k] **z**d + z**d**f **Inv [k] ** j**a + 
z**d**f **Inv [k] **z**d, -j**b - j * *a* * f * * Inv [k] - j**a**g**e + 
j**b**j**b - k**c**h**b - z**d**f **Inv [k] - z**d**g**e + z**e**j**b + 
j**a**f**Inv[k] ** j**b + j**a**f **Inv [k] **z**e + z**d**f **Inv [k] ** j**b + 
z**d**f **Inv [k] **z**e, -a**f**c + b**z**Inv[i] + a**g**d**g**Inv [ i] + 
g**g**g**2**lnv [i] - b** j**a**g**Inv [i] - b* * j * *b* * z * * Inv [ i ] + 
a**g**Inv [i] **h**a**g**Inv [i] + a* *g* * Inv [ i ] * *h* *b* * z * * Inv [ i ] + 
b**z**Inv [i] **h**a**g**Inv [i] + b* * z * * Inv [ i ] * *h* *b* * z * * Inv [ i ] , 
-d**f**c - d**g**Inv[i] + d* *g* *d* *g* * Inv [ i ] + d**g**e**z**Inv [ i ] - 
e** j **a**g**Inv [i] - e** j **b**z**Inv [ i ] - Inv [ i ] * *h* *a* *g* * Inv [ i ] - 
Inv [i] **h**b**z**Inv [i] + d* *g* * Inv [ i ] * *h* *a* *g* * Inv [ i ] + 
d**g**Inv [i] **h**b**z**Inv [i] + e**z**Inv [i] **h**a**g**Inv [i] + 
e**z**Inv [i] **h**b**z**Inv [i] , 

Inv[k]**z - c**h**a**Inv[d] - Inv [k] ** j**a**g - Inv [k] **z**d**g + 
Inv[k] ** j**a**f**Inv[k] **z + Inv [k] ** j**b** j**a**Inv [d] + 
Inv [k] **z**d**f **Inv [k] **z + Inv [k] **z**e** j**a**Inv [d] + 
Inv [k] ** j**a**f **Inv [k] ** j**a**Inv [d] + 
Inv[k] **z**d**f**Inv[k] ** j**a**Inv [d] , 

-Inv[k]**z + c**h**a**Inv [d] - c**h**b**Inv [e] - Inv [k] ** j **b**Inv [e] - 

Inv[k] ** j**a**f**Inv[k] **Inv[e] - Inv [k] ** j **b** j **a**Inv [d] + 

Inv[k] ** j**b** j**b**Inv[e] - Inv [k] **z**d**f **Inv [k] **Inv [e] - 

Inv[k] **z**e** j**a**Inv[d] + Inv [k] **z**e** j**b**Inv [e] - 

Inv[k]**j**a**f**Inv[k]**j**a**Inv[d] + 

Inv[k]**j**a**f**Inv[k]**j**b**Inv[e] 

- Inv[k] **z**d**f**Inv[k] ** j**a**Inv[d] + 

Inv[k]**z**d**f**Inv[k]**j**b**Inv[e], 

''Inv [i] **Inv [c] - Inv[b]**a**f - j**a**g**Inv [i] **Inv [c] - 
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In[3]:= j**b**z**Inv[i] **Inv[c] + z**Inv [i] **h**a**g**Inv [i] **Inv [c] + 

z**Inv [i] **h**b**z**Inv [i] **Inv [c] + Inv [b] **a**g**d**g**Inv [i] **Inv [c] + 

Inv[b] **a**g**e**z**Inv[i] **Inv[c] + 
Inv[b] **a**g**Inv[i] **h**a**g**Inv[i] **Inv[c] + 
Inv [b] **a**g**Inv [i] * *h* *b* * z * * Inv [i] **Inv [c] , 
~z**Inv [i] **Inv [c] + Inv[b]**a**f - Inv[e]**d**f - 

Inv [e] **d**g**Inv [i] **Inv [c] - Inv [b] **a**g**d**g**Inv [i] **Inv [c] - 

Inv [b] **a**g**e**z**Inv [i] **Inv [c] + Inv [e] **d**g**d**g**Inv [i] **Inv [c] + 

Inv [e] **d**g**e**z**Inv [i] **Inv [c] - 

Inv[e]**Inv[i]**h**a**g**Inv[i]**Inv[c] - 

Inv [e] **Inv [i] **h**b**z**Inv [i] **Inv [c] - 

Inv [b] **a**g**Inv [i] * *h* *a* *g* * Inv [ i ] **Inv [c] - 

Inv[b]**a**g**Inv[i]**h**b**z**Inv[i]**Inv[c] + 

Inv[e]**d**g**Inv[i]**h**a**g**Inv[i]**Inv[c] + 

Inv [e] **d**g**Inv [i] **h**b**z**Inv [i] **Inv [c] , 

-d - Inv[i]**h**a + d**g**Inv [i] **Inv [c] **Inv [k] ** j **a + 

d**g**Inv [i] **Inv [c] **Inv [k] **z**d + e**z**Inv [i] **Inv [c] **Inv [k] ** j **a + 

e**z**Inv [i] **Inv [c] **Inv [k] **z**d + 

Inv [i] **h**a**g**Inv [i] **Inv [c] **Inv [k] ** j**a + 

Inv [i] **h**a**g**Inv [i] **Inv [c] **Inv [k] **z**d + 

Inv [i] **h**b**z**Inv [i] **Inv [c] **Inv [k] ** j**a + 

Inv [i] **h**b**z**Inv [i] **Inv [c] **Inv [k] **z**d, 

-e - Inv[i]**h**b - d**g**Inv [i] **Inv [c] **Inv [k] - 

e**z**Inv [i] **Inv [c] **Inv [k] + d**g**Inv [i] **Inv [c] **Inv [k] ** j**b + 

d**g**Inv [i] **Inv [c] **Inv [k] **z**e + e**z**Inv [i] **Inv [c] **Inv [k] ** j **b + 

e**z**Inv [i] **Inv [c] **Inv [k]**z**e - 

Inv [i] **h**a**g**Inv [i] **Inv [c] **Inv [k] - 

Inv [i] **h**b**z**Inv [i] **Inv [c] **Inv [k] + 

Inv[i] **h**a**g**Inv [i] **Inv[c]**Inv[k]**j**b + 

Inv [i] **h**a**g**Inv [i] **Inv[c]**Inv[k]**z**e + 

Inv[i] **h**b**z**Inv [i] **Inv[c]**Inv[k]**j**b + 

Inv[i]**h**b**z**Inv[i]**Inv[c]**Inv[k]**z**e, 

c**h**b + Inv[k]**j**b + Inv[k]**z**e - c**h**a**Inv [d] **e + 

Inv[k] ** j**a**f**Inv[k] - Inv [k] * * j * *b* * j * *b + Inv [ k ] * * z * *d* * f * * Inv [ k ] - 

Inv[k] **z**e** j**b - Inv [k] ** j**a**f **Inv [k] ** j**b + 

Inv[k] ** j**b** j**a**Inv[d] **e - Inv [k] * * z * *d* *f * * Inv [k] * * j * *b + 

Inv[k] **z**e** j**a**Inv[d] **e + Inv [k] * * j **a**f **Inv [k] ** j**a**Inv [d] **e + 

Inv[k]**z**d**f**Inv[k]**j**a**Inv[d]**e, 

c**h**a**g**d - c* *h* *b* * j * *a - Inv [k] * * j * *b* * j**a - Inv [k] **z**d**g**d - 

c**h**a**f**Inv[k] ** j**a - c**h**a**f **Inv [k] **z**d + 

Inv[k] ** j**a**f**c**h**a - Inv [k] * * j * *a* *f * * Inv [k] * * j * *a - 

Inv[k] j**a**f**Inv[k] **z**d + Inv [k] * * j * *a* *g* *d* *g* *d - 

Inv [k] ** j**a**g**e** j**a - Inv [k] ** j**b** j**a**g**d + 

Inv [k] ** j**b** j**b** j**a + Inv [k] **z**d**f **c**h**a + 

Inv [k] **z**d**g**d**g**d - Inv [k] **z**d**g**e** j**a - 

Inv [k] **z**e** j**a**g**d + Inv [k] **z**e** j**b** j**a - 

Inv[k] **j**a**g**d**f**Inv[k] **j**a - Inv [k] ** j **a**g**d**f **Inv [k] **z**d + 
Inv[k] ** j**b** j**a**f**Inv[k] ** j**a + Inv [k] ** j **b** j **a**f **Inv [ k ] **z**d - 
Inv[k] **z**d**g**d**f**Inv[k] ** j**a - Inv [k] **z**d**g**d**f **Inv [ k ] **z**d + 
Inv[k] **z**e** j**a**f **Inv [k] **j**a + Inv[k] **z**e** j **a**f **Inv [k] **z**d. 
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In [4]: = 



c**h**b 
c**h**b* 
2*Inv[k] 
c* *h* *a* 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
-c**h 
Inv [k 
c**h* 
c**h* 
Inv [k 
Inv [k 
2*Inv 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 
Inv [k 



+ Inv[k]**j**b + c**h**a**f **Inv [k] + c**h**a**g**e - 
*j**b + Inv[k] ** j**a**f**Inv[k] + Inv [k] * * j * *a* *g* *e - 
** _ Inv [k] **z**e** j**b - c* *h* *a* * f * * Inv [k] * * j * *b - 

*f**Inv[k] **z**e + Inv [k] * * j * *a* * f * *c* *h* *b - 
j**a**f**Inv[k] ** j**b - Inv[k] ** j**a**f**Inv[k] **z**e + 
j**a**g**d**f **Inv [k] + Inv [k] ** j**a**g**d**g**e - 
j**a**g**e** j**b - Inv [k] ** j**b** j**a**f **Inv [k] - 
j**b** j**a**g**e + Inv [k] ** j**b** j**b** j**b + 
z**cj**f**c**h**b + Inv[k] **z**d**g**d**f**Inv[k] + 
z**d**g**d**g**e - Inv [k] **z**d**g**e** j**b - 
2**e** j**a**f **Inv [k] - Inv [k] **z**e** j**a**g**e + 
z**e** j**]-,** j**b - Inv[k] **j**a**g**d**f**Inv[k] **j**b - 
j**a**g**d**f**Inv[k] **z**e + Inv [k] ** j**b** j**a**f **Inv [k] ** 

j**a**f**Inv[k] **z**e - Inv [k] **z**d**g**d**f **Inv [ k ] * * j ^ 



**Inv [k] **-7**c 
j**a**f**Inv[k] **z**e - 
:**d**g**d**f **Inv [k] **z**e + 
,**o**^**=**f**inv[k] ** j**b + Inv[k]^ 



^b + 

^b - 



*c - 



Inv[i] + Inv[k]^ 



'f**c + 



*f**Inv [k] ^ 



'**d**g**Inv [i] + c* *h* *a* *g* *d* *g* * Inv [ i ] + 



►*Inv [i] 
►*Inv [i] 



* *a* *g* * Inv [ i ] 
Inv [k] **j**a**g**d**f**c 



j * * jj^^-y [ ] + I n V [ k ] * * j * * b * * j * * a 

j**l3**2**Inv[i] - Inv [k] **z**d**g**d 
* * z * *d* *g* *d* *g* * Inv [ i ] - Inv [k] **z**d**g* 
z**e** j**a**f **(, + lnv[k]* 
**a**f**c**h**a**g**Inv[i 



**c 
**c 



^Inv[i] + 



1 

j**a**g**d**g**d**g**Inv[i 
j**a**g**g** j**a**g** jn-j^ [-^ 

j j j **a**g* *Inv [ i 
z**cj**f **c**h**a**g**Inv [i 
2**cj**g**d**g**d**g**Inv [i 
2**(;j**g**e** j **a**g**Inv [ i 
2**e* * j * *a* *g**d**g**Inv [ i 
z**e** j**b** j**a**g**Inv [ i 



j**a**f **c' 



+ Inv [ k ] ** j **a**f **c* *h* *b^ 
+ Inv[k] ** j**a**g**d**g**e* 



►Inv [ 
►Inv [ 



i] + 
i] - 



Inv[k]**j**a**g**e**j 
Inv[k] ** j**b** j**a 



-k r~r-k -k c^-k -k -k -k 



z**Inv [i] 
Inv [d 



i] + 

i] + 

i] + 

i] - 

i] - 

i] + 

+ Inv[k] **z**e** j**b** j**b**z**Inv[i] } 



Inv[k] ** j**b** j**b** j**b**z**Inv[d 
Inv[k] **z**d**f **c**h**b**z**Inv[d 
Inv[k]**z**d**g**d**g**e**z**Inv[d 
Inv[k]**z**d**g**e**j**b**z**Inv[d 



Inv[k]**z**e**j**a 



Inv [d 



(* Here we set the order to be strictly graded lex *) 

SetKnowns[a,Inv[a],b,Inv[b],c,Inv[c],d,Inv[d],e, Inv[e], f, Inv[f], g, Inv[g], h, Inv[h],i, Inv[i],j, Inv[j], k, Inv[k] , z, 
y,x,t,u,v,w]; 

SetUnknowns[{}]; 

(* Ask for a small basis retaining the relations we like *) 
SmallBasis[ allpolys, hopepolys, 3 ] 



Output: 

Out [4]= {— Inv [a] **Inv[h] - Inv[d] **Inv[i] + Inv[a] **b**j** Inv [h] - 



Inv [a] **b**z**Inv[i] — Inv[d] **« 
Inv [d] **e**z**Inv[i], 
■Inv[f] **Inv[a] -Inv[k] **Inv[b] + 



Inv[f ] 



*d * * Inv [a] - Inv [f ] 



'j * * Inv[h] + 



Inv[b] - 



Inv[k] * *z * *d * * Inv[a] + Inv[k] * *z * *e * * Inv[b] } 
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4 Appendix 4 - Confirm Our Relations Imply The Result 



In [5] := Input : 

SetNonCommutative [a,b,c,d,e,f,g,h,i,w,x,y,z,j,u,t,k,v, 
Inv [ a] , Inv [b] , Inv [ c] , Inv [d] , Inv [e] , Inv [ f ] , Inv [g] , Inv [h] , 
Inv[i] , Inv[ j] , Inv[k] ]; 

(* Here we create the relations which we wish to imply *) 
first = {{a,t,b},{u,c,v},{d,w, e}}; 

second = {{ x,f,g}, {h,y,i},{j,k,z}}; 

oneway =MatMult[first,second] - IdentityMatrix[3]; 

otherway = 

MatMult[second,first]- IdentityMatrix[3]; 

start = FIatten[{ oneway, otherway }]; 

SetKnowns[a,Inv[a],b,Inv[b],c,Inv[c],d,Inv[d],e, Inv[e], f, Inv[f], g, Inv[g], h, Inv[h],i, Inv[i],j, Inv[j], k, Inv[k] , z, 
y,x,t,u,v,w]; 

SetUnknowns[{}] ; 

(* hopepolys = *) 

(* hopepolys are defined as in Appendix 3 *) 

(* Here are the relations found through the Small Basis algorithm *) 

newrels = {-Inv[a]**Inv[h] - Inv[d]**Inv[i] + Inv[a]**b**j**Inv[h] - 
Inv[a]**b**z**Inv[i] - Inv[d]**e**j**Inv[h] + Inv[d]**e**z**Inv[i], 
-Inv[f]**Inv[a] - Inv[k]**Inv[b] + Inv[f]**g**d**Inv[a] - 
Inv[f]**g**e**Inv[b] - Inv[k]**z**d**Inv[a] + Inv[k]**z**e**Inv[b]} 

hopepolys = Join[ hopepolys, newrels ]; 

(* We will have our result if the starting relations are elements of 
the ideal generated by the above relations i.e. the GB created with 
the above relations reduces the starting relations to *) 

(* Note that a 5 iteration GB failed to reduce the starting relations *) 
hopeGB = NCMakeGB[ hopepolys , 7 ]; 

hoperules = PolyToRule[ hopcGB ]: 

(* Use the Groebner basis relations to reduce the original relations *) 
Reduction! start , hoperules ] 

Output: 

Out[5]= {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} 
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